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Abstract 

We compute the contribution of the vacuum Virasoro representation to the genus-two 
partition function of an arbitrary CFT with central charge c > 1. This is the perturbative 
pure gravity partition function in three dimensions. We employ a sewing construction, in 
which the partition function is expressed as a sum of sphere four-point functions of Virasoro 
vacuum descendants. For this purpose, we develop techniques to efficiently compute correla¬ 
tion functions of holomorphic operators, which by crossing symmetry are determined exactly 
by a finite number of OPE coefficients; this is an analytic implementation of the conformal 
bootstrap. Expanding the results in 1/c, corresponding to the semiclassical bulk gravity 
expansion, we find that—unlike at genus one—the result does not truncate at finite loop 
order. Our results also allow us to extend earlier work on multiple-interval Renyi entropies 
and on the partition function in the separating degeneration limit. 
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1 Introduction 

Three-dimensional gravity has proven to be a frnitfnl testing gronnd for onr ideas abont 
holography. Perhaps the most interesting qnestion is whether pnre general relativity—a 
theory with only metric degrees of freedom—with a negative cosmological constant exists as 
a qnantum theory in its own right. If this were the case, then one should be able to find 
its holographic dual for a given value of GN/RAds, the Newton constant in AdS units. This 
appears to be an extremely difficult problem (see e.g. Si)- However, general relativity 
exists as a sub-sector of any theory of gravity in three dimensions. From the boundary 
point of view, it captures the dynamics of the Virasoro sector of any two-dimensional CFT 
with central charge c = 3i?Ads/2G'Ar- This semi-microscopic interpretation is unavailable in 
higher-dimensional AdS/CFT, where the stress tensor does not generate a closed symmetry 
algebra. 

This perspective lends a universality to AdS 3 /CFT 2 that underlies, for example, the 
matching of the asymptotic symmetry group of anti-de Sitter space to the Virasoro algebra 
and the matching of BTZ black hole entropy to the Cardy growth of states [^. These are 
features of any theory of three-dimensional gravity, and of any dual CFT. Recent work has 
revealed an even richer set of properties of two-dimensional CFTs that admit a large-c limit 
and are dual to weakly coupled bulk theories of gravity. These relate to aspects of such 
theories’ spectra and thermodynamics [6-10 , entanglement 11-19 , Virasoro blocks 20-24 


modular geometry 25 , and chaotic response to perturbations 26 , among others. 


Despite this fascinating progress, much remains to be understood about basic conse¬ 
quences of Virasoro symmetry. To this end, in this paper we will focus on the computation 
of the partition function of three-dimensional gravity in a universe whose boundary is a 
Riemann surface of genus g. Schematically, this should be given by a bulk path integral 
over geometries Aig which asymptote to 


Z^r.A^g) = / Vg . (1.1) 

JdMg=i:g 

This partition function is a function of the conformal structure moduli of the Riemann surface 
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Eg, denoted ^Ig. These partition functions contain vital information about the theory: for 
instance, one can recover the correlation functions of a given CFT from its higher-genus 

Thus, by tuning the moduli Qg, one could in 


partition functions by pinching handles 27 


principle recover the correlation functions of the boundary CFT. 


Equation (1.1) is in general an extremely difficult object to compute. Moreover, it is not 
■‘universal” in the sense described earlier. In particular, in (1.1) we have written the bulk 


path integral only over metric degrees of freedom; in more complicated theories of gravity 
more degrees of freedom should be included. The partition function Zgj.a.v{flg) written above 
is that of the CFT dual to pure gravity at a given value of Newton’s constant, if it exists. In 


this paper we will not be interested in the full partition function (1.1), but rather in an object 
which is both easier to compute and universal: we will study the contribution to .^grav(f^g) 
from a single saddle-point geometry Aig, including perturbative quantum corrections. This 
restricted partition function maps to the contribution of the Virasoro sector to the CFT 
partition function on 

This is easiest to see at genus one. In the semiclassical regime Gat <C i?Ads, fhe path 


integral (1.1) can be recast as a sum over saddle points of the Einstein action with solid torus 


topology, along with perturbative corrections. The simplest such saddle point is thermal 
AdSa, the Euclidean geometry found by taking empty AdSa and periodically identifying in 
Euclidean time, which contributes to the partition function as 


ZTAds{r,T) = \q 


_ 1 ^ 1 -c/12 


n 

n=2 


12 ’ 


q := e 


^2'KiT 


( 1 . 2 ) 


In this expression, we have included not only the classical action of thermal AdS (the factor 
of but also all of the perturbative quantum corrections which come from loops of 

gravitons in thermal AdS. With certain reasonable assumptions, all other saddle points are 
simply S'L(2,Z) modular transformations of thermal AdS, and the sum over geometries is 


a sum over S'L(2,Z) transformations of (1.2). A direct calculation of .^xAdS does not yield 
a result consistent with an interpretation as a trace of the Hilbert space of a CFT iHIQ 


Nevertheless, (1.2) does have a natural interpretation as the Virasoro vacuum character of 
any CFT with central charge c > 1 and an S'L(2,M) x S'L(2, M)-invariant ground state. We 


note that this object is not modular invariant, which reflects the fact that in (1.2) we have 


focused on only one saddle out of the SL{2, Z) family. In the language of Riemann surfaces. 


(1.2) is a function not of the conformal structure of the boundary torus, but rather of the 


Teichmiiller parameter r. 


In the present paper, our goal is to compute the analog of (1.2) at higher genus. Any the- 


^However, in the quantum regime Gn ^ /?AdS) it was argued in that at specific minimal-model values 
of c, not only can the sum be performed, but it agrees with the minimal-model partition functions. 
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ory of AdSa gravity contains solutions which are handlebodies—solid genus-^f geometries— 
which have the Riemann surface Eg as their conformal boundaries. These solutions are quo¬ 
tients of Euclidean AdSa, much like thermal AdS. The contribution of a given handlebody to 
the path integral—including graviton loop corrections—has a universal CFT interpretation 
for any value of c, as the contribution of the states in the vacuum representation to the 
CFT partition function on Tig. We call this quantity 2’vac0 Z^ac is a function not just of 
the conformal structure of Tg, but rather of the Teichmiiller parameters that parametrize 
the universal cover of the moduli space. In other words, to compute Zvac we must specify a 
marking of the Riemann surface Tg that Exes a choice of contractible and non-contractible 
cycles of the handlebody (A- and B-cycles, respectively). Thus Zyac is not invariant under 
the modular group; a modular-invariant partition function could be obtained only, for ex¬ 
ample, by summing over bulk saddles that describe the different ways a boundary Tg can be 
“hlled in” by bulk geometries. One of our goals in this paper will be to give a direct CFT 
computation of ^vac, which can then be interpreted gravitationally. 

There has been a recent resurgence of interest in higher-genus partition functions of two- 
dimensional CFTs. This interest is partly motivated by the study of entanglement entropies 
(EEs). The computation of EEs via the replica trick involves evaluating entanglement Renyi 
entropies (EREs), which in turn are equal to certain higher-genus partition functions. A 
particularly interesting line of research uses calculations of EREs to test the Ryu-Takayanagi 
(RT) classical EE formula and understand the quantum corrections to it 

The partition functions relevant for EREs have been computed in holographic CFTs in 
two ways: in gravity, by explicitly Ending the relevant saddles and evaluating their classical 
actions and one-loop determinants, and in held theory, by computing twist-operator corre¬ 
lators in certain cyclic orbifold CFTs and then expanding the results in powers of 1/c. In 
every case where the computation has been carried out on both sides, agreement has been 
found. This is a check of our basic understanding of AdS 3 /CFT 2 duality. Further, in many 
cases the computation using one technique gives results that are not practically computable 
using the other, thereby giving new information about both three-dimensional gravity and 
large-c CFTs. For example, by expanding the results of the twist-operator computation to 
higher orders in 1/c, one determines higher-loop quantum corrections on the gravity side 
that would be exceedingly difficult to obtain by direct computation. These results hint at a 
surprising novel structure, which we describe below. 


11,13-15 


^This was called Zfake in 
written as Zfake = Zsaddle- 


28 , and the correspondence with the bulk saddle point partition function was 
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Figure 1: A depiction of the sewing construction as applied to ^vac, the contribution of the 
Virasoro vacuum representation to a genus-two CFT partition function. The coordinates pi 
and p 2 represent the widths of the two handles in a Schottky uniformization of the Riemann 
surface. The handles are replaced by a sum over pairwise operator insertions, where we 
include all Virasoro descendants of the identity, O G "Hvac- This recasts Zvac as a sum of 
sphere four-point functions, weighted by powers of pi and p 2 - The operators Ot and Oj have 
holomorphic conformal weights hi and hj, respectively. A detailed description of the sewing 
construction is presented in section]^ (see equations 4.12 and 4.13). 


1.1 Summary of results 


In this paper we will directly compute Zvac at genus two, for arbitrary values of c > 1, using 
CFT techniques. We will use a sewing construction, represented schematically in £gure[T] We 
start with a Riemann surface S that has been constructed by Schottky uniformization and 
replace the handles of S with a sum over states propagating along these handles. The result 
is a weighted sum over four-point functions of local operators on the sphere. If we were 
to include all possible operators in this sum, we would obtain the full, modular-invariant 
partition function, as a function of the pinching parameters pi and p 2 , that describe the 
widths of the handles, and a third modulus x, which is the cross-ratio of the four-point 
functions on the sphere. The universal contribution Z^ac is computed by summing only over 
operators in the Virasoro vacuum block. The four-point functions of these operators are 
determined completely by conformal Ward identities. Thus Zvac is in principle completely 
determined in terms of the central charge. We will compute the answer perturbatively 
in Pi and p 2 but exactly in x. We will mostly assume that the CFT has no extended 


symmetry algebra beyond two copies of Virasoro, although as we will see in section 5.1.1 it 
is straightforward to extend our results to higher-spin theories. 


Conformal bootstrap methods play an important role in our computation of Zv 


since 


our computation requires us to sum over all four-point functions of Virasoro descendants 
on the sphere. These correlation functions—and indeed the correlation function of any 
family of chiral operators—can be efficiently computed using a holomorphic version of the 
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conformal bootstrap. The essential idea is that these correlation functions can be regarded 
as a meromorphic functions on A^o ,45 the moduli space of four marked points on a sphere, 
with poles only when the operators coincide. A meromorphic function on a compact space is 
determined entirely by its polar behaviour. For chiral operators of hnite conformal dimension, 
this polar part is determined by a hnite number of three-point function coefficients. The 
result is an exact expression for the correlation function in terms of a finite number of 
three-point function coefficients. This should be contrasted with the usual approach, where 
a four-point function involves a sum over an inhnite number of intermediate states, so is 
written in terms of an inhnite sum of OPE coefficients. Similar ideas have been advanced 


m 


30-32 . When the chiral operators are Virasoro descendants of the identity, we show 


using free bosons that all connected n-point functions have polynomial dependence on c. 
This implies that, when expressed in terms of c, bulk scattering of graviton states in AdSa 
is purely classical, in analogy to the one-loop exactness of the torus partition function. 

Our result for will hold for a general Riemann surfaces S, but for certain values of 
the moduli—those corresponding to the so-called replica surface—our results can be used 
to compute genus-two EREs. We mainly consider the case of two disjoint intervals in the 
vacuum of a GET; the replica manifold has genus two when n = 3, and is denoted <^ 2 , 3 - Our 
results extend previous results in 16,33 , which were obtained from the twist-held four-point 
function. Those works employ a short-interval expansion in which the conformally invariant 
cross-ratio, which we call y, is taken to be small. The sewing technique is well-suited to 
computation to higher orders in y, 33 worked through 0{y^), and we extend this to 0{y^^). 
In fact, the authors of 33 found a quite remarkable result: their 0{y^) term in logZ^ac 
exhibits a two-loop truncation in the expansion in 1/c at large c. 

To understand why that result is interesting, let us consider the bulk AdSs interpretation 
of our results. Our computation of Zvac is not limited to large c; it is a truly quantum result 
for the saddle-point partition function for a genus-two handlebody of three-dimensional pure 
gravity, applicable even when Gat/RacIS is of order one. Expanding our result at large c 
is equivalent to making the semiclassical approximation in the bulk. More precisely, the 
expansion of the “vacuum free energy,” 


Tvac • log Zvac ; (i- • ) 

at small Newton constant (large c) is the loop expansion of three-dimensional AdS gravity: 

00 

R.ac = (1.4) 

i=0 

where denotes an Gloop contribution. In the bulk, no computations have been done 
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beyond one-loop order. At one loop, there is a closed-form expression for the graviton 
handlebody determinant 28,34 . Our result for Fvacp is a computation of this determinant 


in a new regime of moduli space, not described by previous computations 15,28 


What about at higher loops? At genus one, the expansion ( |1.2[ ) truncates at one-loop 
order: higher-loop contributions only renormalize the value of c [^. It is natural to ask 
whether the higher-genus partition function obeys an analogous truncation. Indeed, the 
results of 33 imply that that Fiy 2 {^ 2 ,z) vanishes through perhaps suggesting that 

the partition function at genus two truncates at two loops. One motivation for this paper was 
to investigate whether this truncation really occurs for the full partition function Our 

conclusion is that the truncation does not occur, and that the cancellation observed in is 
an artifact of the small-y expansion. Indeed, we will show that on the replica manifold <^^ 2,3 
there are nonzero contributions to the free energy at all orders in the 1/c expansion. These 
hrst appear at in the short-interval expansion, explaining why these corrections were 

not found in (^ . 

More generally, we will show that the genus-two partition function Zvac of pure three- 
dimensional gravity does not truncate at any order in 1/c. The same is true of pure higher 
spin gravity. Explicit contributions to the all-loop terms Evac;£ are given in section 5.1 To 


our knowledge, these are the first all-loop results beyond genus one for a Riemann surface 
with three independent moduli. We show that in the regime of small pi and p 2 i fhe only 
point in the moduli space at which the loop expansion (1.4) truncates is the separating 
degeneration point, at which S degenerates into the union of two tori. 

This paper is organized as follows. In section we recall the relationship between Renyi 
entropies and higher-genus partition functions, and review the sewing construction of higher 
genus partition functions as a weighted sum over sphere-correlation functions. In section 
we describe techniques to compute these correlation functions, including an analytic version 
of the conformal bootstrap. In section we apply these techniques to compute .^vac, the 
contribution to the partition function from the Virasoro descendants, at genus two. We 
discuss the large central charge limit of this result in section which allows us to understand 
the nature of quantum corrections to the higher-genus partition function of three-dimensional 
gravity, as well as applications to Renyi entropies, before concluding in section]^ Appendices 


^ On general grounds, such a truncation might seem to conflict with the pole structure of CFT correlation 
functions, regarded as analytic functions of c. Let us make the argument at genus two for concreteness. 
In the sewing construction, a genus-two partition function is written as a sum over four-point functions. 
The statement of truncation becomes the statement that at each order in the sewing expansion, the total 
contribution from all four-point functions truncates at order 1/c. As argued by Zamolodchikov 35 , the 


conformal block decomposition of a given four-point function contains poles at minimal model values of c 
where the exchanged operators become null. Unless these poles cancel against the poles in the other four- 
point functions contributing at a given order in the sewing expansion, the partition function will not truncate 
in a 1/c expansion. This sort of cancellation at every order in the sewing expansion seems highly unlikely. 
Indeed, our computations bear out this conclusion. 
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contain details relevant to the sewing construction. 


2 Review 


In this section we will review some relevant background material, and explain the methodol¬ 
ogy and philosophy behind our computations. In subsection |2.1 , we briefly review previous 
work on Renyi entropies in the vacuum of 2D CFTs. Using these Renyi entropies as a guide, 
we explain how pure 3D quantum gravity naturally computes the universal contribution of 
the Virasoro identity block to CFT partition functions on generic Riemann surfaces. Then, 
in subsection 2^, we will explain the sewing construction, which we will apply in section]^ 
to the computation of higher-genus partition functions. 


2.1 Renyi entropies and higher-genus partition functions 

Two-dimensional CFTs provide perhaps the simplest arena in which to investigate entan¬ 
glement entropies (EEs) in held theories. In this subsection, we will briehy review some 
calculations of these quantities, with particular attention to their dependence on the central 
charge c of the theory. 


2.1.1 General CFTs 

The simplest quantity one can consider in this context is the vacuum EE of a single interval 
[m, n] on the line. The corresponding Renyi entropy is given in terms of the partition function 
on the surfac^.^i^„, which is the plane branched n times over the interval 36 

1 




n — 1 




l,n ) 


( 2 . 1 ) 


This partition function is in turn related to the two-point function on the plane of twist 
operators in the orbifold theory C”/Z„ (where C is the original CFT) (^: 




( 2 . 2 ) 


It will be convenient to work in terms of the free energy, which we dehne on any surface X 


as 


F{X) := -lnZ(X). 


(2.3) 


^In general, we will denote the plane branched n times over a set of N intervals by ^N,n', this surface has 
genus {N — l)(n — 1). 

®The partition function on a genus-zero surface is defined, for a given theory, up to a multiplicative 
constant independent of the metric. We are choosing that constant so that Z{C) = 1; otherwise the argument 


of the logarithm in (2.1) would be Z(.^i^„)/Z(C)". 











The free energy is proportional to c and otherwise independent of the theory: the 

surface has genus zero, so the free energy is given entirely by a Liouville action multiplied 
by c; alternatively, in the twist-operator language, the two-point function depends only on 


their dimension, which is proportional to c. The result is 36,37 




n 


(2.4) 


where e is an ultraviolet-cutoff length scale; its presence reflects the divergence in the par¬ 
tition function due to the presence of conical singularities in Mi^n- This gives rise to the 


well-known formula for the EE 36,37, 


^([M,n]) = = 3 In 


V — u 


(2.5) 


The above result is easily generalized to the case of a single interval on a circle at zero 

In either case, the branched cover 


temperature or on the line at hnite temperature 37 


surface continues to have genus zero, and therefore the EREs and EE continue to be propor¬ 
tional to c. The simplest cases where higher-genus partition functions appear are a single 
interval on the circle at hnite temperature and two intervals on the line at zero temperature; 
the corresponding branched-cover surfaces have genera n and n — 1, respectively. This im¬ 
plies that the ERE will depend on the full operator content of the GET, not just its central 
charge 38 . In rest of this subsection we will focus on the two-interval case, which is the 


best-studied one. 

For two intervals [Mi,ni] U [u 2 -iV^^ it is convenient to work with the mutual (Renyi) 
information, which is ultraviolet-hnite, hence conformally invariant and dependent only on 
the cross-ratio y of the four endpoints 38 


/'">(</) := S<">(|'t.l, !.,|) + S'">(|«2. !>2|) - S<">(|'«1. !>l| U [«2, t>.2]) , V ■= _ “jj ' 

"Le) 


We have 


- F{^ 2 ,n) -l-subtractions . 


1 — n 


(2.7) 


The subtractions, given by the EREs of the individual intervals, soak up the divergences in 
F{^ 2 ,n), leaving an unambiguous hnite value for The partition function on ^ 2 ,n can 

be expressed as a four-point function of twist operators in the orbifold theory: 


Z{^2,n) ^ {(^{ui)r^[Vl)(T{u2)a{v2))c’'j7^ 


( 2 . 8 ) 


®In the literature, this cross-ratio is often denoted x; however, we will use x for a different cross-ratio in 
what follows. 
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Figure 2: The u-sheeted replica surface which is the branched covering surface of the 
plane with two intervals and has genus n — 1. On each sheet, there is a cycle separating the 
two intervals called the A-cycle, and another cycle encircling the two points Vi and U 2 -, called 
the -B-cycle. There are n — 1 independent cycles of each type. 

The surface has genus n — 1, so the partition function depends on the full operator 
content of the theory and not just its central charge. However, it contains a universal 
contribution that only depends on c. To dehne this part, it is useful to hrst set up some 
notation regarding the topology of the surface ^ 2 ,n- 

A useful basis of cycles on ^ 2 ,n can be described as follows. On each sheet, there is a 
cycle that separates the two intervals. We will call these A-cycles. The sum of all n of them 
is trivial, so there are n — 1 independent ones. There are also cycles which encircle the points 
ni,M 2 , crossing each cut once, which we call B-cycles; again, there are n — 1 independent 
ones. (See £gure|^) The A-cycles intersect the B-cycles but not themselves, and vice versa. 
(Linear combinations A,, Bj can be constructed with intersection numbers Aj ■ Bj = Stj, but 
this will not be necessary for our purposes.) It is also useful to visualize the surface ^ 2 ,n 
as two spheres connected by n tubes. This can be related to the branched cover by cutting 
each sheet along a small ellipse surrounding the interval [ui, ui] and another one surrounding 
the interval [u 2 ,V 2 \- Each interval then becomes a sphere with n holes, while each sheet 
becomes a tube connecting one sphere to the other. Each A-cycle wraps a tube, while each 
B-cycle runs along one tube and back along another. (See £gure[^) ^% 2 ,n enjoys a “replica 
symmetry”, which cyclically permutes the sheets, and hence also the tubes. 

The universal part of Z{M 2 ^n) to which we alluded above is defined as the contribution 
in which only Virasoro descendants of the vacuum appear on the A-cycles. In other words, 
if for any circle C we dehne Pvac(C) as the projection operator onto the conformal family of 
the vacuum of the Hilbert space of C on C, then we dehne the vacuum partition function as 
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A 



Figure 3: An alternate depiction of the surface ^ 2 ,n in figure]^ ^ 2 ,n can be visualized as two 
spheres connected by n tubes. The two spheres, one for each interval, are made by cutting 
small holes around each pair of intervals on all n sheets. The tubes connecting the holes on 
the two spheres represent the sheets. In this picture, the A-cycles wrap the n tubes and the 
5-cycles run through two different tubes. 


the path integral with projectors Pvac(Ai) ■ ■ ■ Pvac(An_i) inserted: 

^vac('^2,n) := (Pvac(Ai) ■ ■ ■ Pvac(A„_i))Z(^2,n) • (2-9) 

With only vacuum descendants on Ai,..., A„_i, the cycle A„ (which is a linear combination 
of the others) is automatically guaranteed to host only such descendants as wellj^ Note that 
the choice of representative of any given cycle A* is unimportant; any representative can be 
mapped to any other by a holomorphic diffeomorphism, which acts on the Hilbert space by 
the Virasoro group, under which conformal families don’t mix. In the orbifold description, 
the vacuum partition function can be written 

•^vac('^2,n) (cr(^l)^(^l)Porb vac(A)(T('U2)<^(^2))c"/Zn ) (2T0) 

where Porb vac is the projector onto states of C"/Z„ composed of descendants of the identity of 
C|^and A is a circle enclosing [ui, ui]. Note that, unlike the full partition function, .^vac('^ 2 ,n) 
is not a modular invariant quantity, due to the distinguished role of the A-cycles. 

As we will see in the next subsection, the vacuum partition function is particularly well- 

^To see this, cut along all n A-cycles, leaving two sphere n-point functions (on the left and right spheres 
of figure]^. For each n-point function, n — 1 of the operators correspond to vacuum descendants. As a 
result, if the nth one is not a vacuum descendant, then the n-point function vanishes, and hence does not 
contribute to the vacuum partition function. 

®This set of states includes more than just Virasoro descendants of the vacuum of C”/Z„. Rather, it 
includes all descendants of the vacuum under the larger algebra consisting of (Z„-symmetric) products of 
Virasoro generators acting on the different copies of C. 
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studied in the context of holographic and other large-c CFTs. 


2.1.2 Large-c CFTs 

We are interested in families of CFTs, such as holographic ones, that admit a large-c limit. 
In such theories, all of these quantities—the free energies, entanglement (Renyi) entropies, 
and mutual (Renyi) informations—admit an expansion in 1/c starting at order c. We thus 
write, for example, 

OO CXD 

, F(^2,n) = • (2.11) 

1=0 £=0 


In a holographic CFT, the parameter 1/c is proportional to the bulk Newton constant. 


2G 


N 


c SRAds ’ 


( 2 . 12 ) 


so the expansion in 1/c is a loop expansion in the bulk (hence the index i). 

From a CFT perspective, the loop corrections {£ > 1) are “cleaner” than the classical 
one (£ = 0), in the following sense. First, F/>i is unambiguous, since the scheme dependence 
of the free energy is due to the Weyl anomaly, which is proportional to the central charge. 
Second, it is hnite even on a singular surface such as ^ 2 ,n, since the Weyl transformation 
that smoothes out those conical singularities shifts the free energy by c times a Liouville 
action. Third, since it is Weyl-invariant, it depends only on the complex structure of £^ 2 ,n, 
hence only on the cross-ratio y, not the positions of the endpoints themselves. Finally, since 


the subtractions present in (2.7) are proportional to c, we simply have 


G(y) = 


1 — n 


Fi{^2,n) 


for £ > 1. 


(2.13) 


These properties will all be useful when we study the loop corrections below. 

The RT formula makes a strikingly simple prediction for the classical part of the mutual 


information 11 


Io{y) = 


0 , 2 /< 1/2 

(l/3)ln(2//(l-2/)) 2 />1/2 


(2.14) 


It is interesting that this formula does not depend on the held content or other specihcs of 
the dual theory. On the other hand, the loop corrections do depend on the held content, 
although they always include certain “universal” terms due to the gravitational sector, as 
we will explain below. 


Signihcant ehort has gone into testing the prediction (2.14) and computing the loop cor 
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rections using the replica trick. Two strategies have been followed to compute the relevant 
free energies. The hrst is to hnd the dominant gravitational saddle point whose conformal 
boundary is the terms in the 1/c expansion of F{^ 2 ,n) are then given by the classi¬ 

cal action, the one-loop determinant of the helds about that background, and so on. The 
second strategy is to compute the four-point function of twist helds (2.8) using CFT tech¬ 
niques such as the conformal-block decomposition. The RT prediction for the classical part 


was successfully conhrmed, modulo some assumptions, by both methods, in 14 and 13 
respectively. 

Consider the calculation of id)('^ 2 ,n), starting with the gravity method. In 14 , two grav¬ 
itational saddles were constructed with conformal boundary ^ 2 ,n- Both are handlebodies; in 
one, which we will call Ha, the A-cycles are contractible, while in the other, Hb, the B-cycles 
are contractible. Ha has a smaller action for y < 1/2 and Hb for y > 1/2. These are the 
only solutions that preserve the replica symmetry of ^% 2 ,n, and are also believed to be the 
only type of solution that exists uniformly for all n. Their actions are analytic functions of 
n; when continued down to n = 1, they reproduce precisely the RT prediction (2.14) for the 
EE. In 13 , an analysis of conformal blocks in the theory at large c—again, imposing 

the replica symmetry—led to the same result. 

An important subtlety regarding these calculations is as follows. For general n and y, it 
is not clear whether the dominant gravitational saddle is always Ha oi Hb, and therefore 
whether their actions indeed give the correct free energy and Renyi entropy]^ However, for 
small y, the tubes are very thin (as we will see in section 5.2 when we discuss the period 
matrix for so the dominant saddle must indeed be the handlebody that hlls them in, 

namely Ha- This is important for our purposes because the calculations we will describe 
from here on will always be done in an expansion in ?/, and therefore we can safely ignore 
this subtlety and assume that Ha is the dominant saddle. 

We now turn to the one-loop correction to the free energy, Fi[M 2 ^n), which as noted 
in (2.13) directly gives the one-loop correction to the mutual Renyi information, I^\y). 
Fi{^ 2 ,n) is proportional to the sum of the logs of the fluctuation determinants of all the 
helds propagating on the relevant gravitational saddle. In any theory of gravity, this includes 
the metric huctuations. Their one-loop determinant on the handlebody Ha was computed 
in an expansion in y to order y^ for all n in 15 , and to order y^^ for n = 1 in 1^. 


As we will now explain, this contribution to the free energy is simply the 0{c^) part of the 
vacuum free energy Z^s_c{^ 2 ,n) dehned in the previous subsection. In fact, more generally, 
consider the partition function obtained from the classical action and loop corrections to 
all orders of perturbative pure gravity on Ha- We will now argue that this quantity is 


®Even if this is not the case, one can argue that these are the relevant saddles to consider for the purposes 
of analytically continuing the ERE down to n = 1 to find the EE. 
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precisely ^vac('^ 2 ,n)- In the genus-one case, this was shown in |^, and we can adopt their 
argument here. In a Hilbert-space interpretation, we can choose to think of the A-cycles as 
dehning a spatial direction and the B-cycles a (Euclidean) time direction. This is convenient 
because the states dehned on the A-cycles are perturbative pure quantum gravity states 
on an AdSs background, since the A-cycles are contractible and the handlebody is locally 
Euclidean AdSs. Since the creation operators for metric fluctuations are, from the GET 
viewpoint, Virasoro generators, these states are thus Virasoro descendants of the vacuum. 
Thus the perturbative pure gravity partition function on Ha is precisely 2’vac(^2,n)- The 
exact correspondence between the perturbative quantum gravity partition function and the 
universal identity block contributions to GET partition functions was articulated and tested 
in 


28 . We will extend that work in section 5.3 


We now return to computation of Renyi entropies specihcally. Having established the 
GET interpretation of ^vac('^ 2 ,n), we can see that reproducing the results of 15,39 using the 


twist-held method requires including only descendants of the vacuum as intermediate states 


in the conformal-block decomposition of the 4-point function (2.8), since the intermediate 


states are precisely those living on the A-cycles. More precisely, one should include states of 
the orbifold theory that are made up of descendants of the vacuum of C; these include 

more than just the descendants of the vacuum of j^Ln- It is easy to see that the term of 
order in is given by descendants at level h. 

These calculations were carried out to order y^ by Ghen et ah in 33 . Expanding their 
result in powers of 1/c, the one-loop (order c°) term matched the bulk metric one-loop 
determinant computed earlier in 15 . Their term started at order so their results 


could access i < 3. In other words, they not only reproduced the one-loop determinant, 
but effectively computed two-loop and three-loop free energies, which would presumably be 
quite challenging from a direct bulk perturbative calculation. The coefficient at each order 
in 1/c and ?/ is a rational function of n. We will not reproduce these rather complicated 
functions here for general n. However, let us note the following pattern in the n-dependence 


observed by 33 : 




c,£(,^ 2 ,n) = {n- l)(n - 2) ■ ■ ■ (n - ^) ( ^ amXn')p^ ) + 0(:X') 


(2.15) 


\ m = 2£+2 


The cXm/in) are functions of n; some of them have zeroes at positive n, but none of these 


zeroes coincide, unlike those shown in (2.15). 


There are some notable features of this formula. First, Fvac; £>i(^ 2 ,n) carries a factor of 
n — 1. The fact that it vanishes at n = 1 can be understood from the fact that the genus- 
zero free energy is given entirely by a Liouville action multiplied by c; it is also necessary. 
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given (2.13), for the mutual Renyi information to have a smooth limit as n —)■ 1. Second, 
Fvac-, £> 2 {^ 2 ,n) Carries an overall factor of n — 2. The fact that it vanishes at n = 2 can be 
understood from the fact that the contribution of the identity family to the genus-one free 
energy is one-loop exact: aside from a classical (order-c) part, it is given by — lnxvac(?/), 
where Xva.c{y) is the character of the identity family, which is independent of c. 

Perhaps surprisingly, the j/® term of TVac; 3 ('^ 2 ,n) computed by Chen et ah carries an 


overall factor of n — 3. (Recall that 33 only computed through 0{y^) in the |/-expansion.( 


If the pattern (2.15) were to hold to all orders in ?/, this would imply a truncation in the loop 


expansion around handlebodies asymptotic to ^ 2 ,n with appropriate cycles contractible. On 
this basis, Chen et al. were led to suggest that the genus-^f free energy might be ^f-loop exact 
for all g, at least for the replica manifolds ^N,n- One might even wonder whether this could 
be true for all genus-gf manifolds. One of the main purposes of this paper is to test this 
intriguing idea. To do this, we will calculate Z^ac on generic genus-two Riemann surfaces, S, 
using a different technique that we describe now; this complementary approach will provide 
a gateway to applications to Renyi entropy and 3D quantum gravity. 


2.2 Vacuum amplitudes from sewing 

In section we will compute Zvac via the sewing construction. We heuristically explain 
this method here with the help of figure the method applies to computation of the full 
partition function Z of C, but can be specialized to computation of Zvac- The basic idea 
is to replace each handle of S by a sum over local operator insertions at its ends. This 
frames the computation of Z as a weighted sum of sphere four-point functions. As stressed 
earlier in this section, computing Zvac as opposed to the full partition function of C means 
that we only allow Virasoro vacuum descendants to propagate along the handles. This 
construction is perturbative in the width of the handles. There are many parameterizations 
of a given surface S; we use the Schottky construction, which forms S as a quotient of 
the Riemann sphere by a discrete subgroup of PSL{2, C), the Mobius group. The genus-two 
Schottky space is parameterized by coordinates {pi,p 2 ,x}. Roughly speaking, these describe 
the width of the two handles and the sphere coordinate of the lone endpoint not fixed by 
conformal symmetry, respectively. The computation of Zvac is then a double power series in 
Pi and p 2 , where the powers are the left-moving conformal weights of the operators inserted 
at the endpoints. 


In order to make eventual contact with Renyi entropies and the work of 28 , we will also 
need to express Zvac in terms of the period matrix of S, denoted D. That is, we need to 
perform the coordinate map {pi,P 2 , a:} t D. This is known in closed form, but is complicated 
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Figure 4: A picture of the sewing approach to computing a genus-two partition function, 
Z. The mechanism was explained in hgure [T] To compute Z rather than one simply 
lets the sum run over all operators in the CFT Hilbert space. 


see, 


, e.g., 28,^). If we define multiplicative periods 


qij := e 




then Qij admits a power series in pi and p 2 of the following form: 

oo n+m 

qii=Pi'^PlP 2 X] c{n,m,\r\) x\ 

—m 
n+ra 

d{n,m,r):j 


qi 2 = X + X ^ 
n,m=l 


P>2 


r=—n—m 


922 = 9ll(Pl P 2 ) ■ 


(2.16) 


(2.17) 


The c(n, m, |r|) and d{n,m,r) = d{m,n,r) are coefficients given in Appendix E of 40 
through m = n = 7. 

Thus, in order to compute Zvac via sewing, we must compute four-point functions of 
operators in the Virasoro identity representation. We turn to this now, by way of the more 
general subject of computing four-point functions of arbitrary holomorphic operators. 


3 Four-point functions and the analytic bootstrap 

In this section, which may be read independently of the rest of the paper, we discuss methods 
for computing correlators in 2D CFTs. We will focus on four-point functions. A standard 
way to compute a four-point function is to do an OPE expansion of pairs of operators. This 
yields a power series in the cross-ratio x of the four points. However, we wish to calculate 
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the correlator at finite values of x. We will describe two methods to do this. The first, 
described in subsection|3.2[ is via direct manipulation of operator modes, and may be familiar 


to CFT practitioners. The second, described in subsection |3.3[ is an analytic realization 
of the conformal bootstrap that applies specifically to correlators of only holomorphic (or 
only anti-holomorphic) operators. The upshot is that the combination of holomorphy with 
fundamental properties of the OPE and crossing symmetry yields a result that is determined 
by a finite number of OPE coefficients. In the case that all four operators have identical 
holomorphic dimensions, the solution of crossing symmetry leads to an especially simple 
algorithm. 

Before turning to those methods, we first describe a simple approach that applies specif¬ 
ically to correlators of descendants of the identity, and use the results to discuss the powers 
of the central charge c that appear in such correlators. 


3.1 Free bosons and powers of c 

As explained in subsection (and illustrated in figure [^, Zvac, the universal part of the 
genus-two partition function, can be constructed from four-point functions on the plane of 
descendants of the identity. In this subsection, we will discuss general properties of such 
correlators. 

The first property to note is that they are independent of the rest of the field content of 
the theory, and depend only on its central charge. This follows from the fact that the identity 
Verma module is closed under fusion. Since in this paper we are particularly interested in 
the powers of c that appear in Zyac, in this subsection we will focus on the question of what 
powers of c can appear in such correlators. We will first use a simple counting argument 
in a free field theory to show that the powers of c are highly constrained. In particular, if 
one thinks of 1/c as a coupling constant, then it appears that these correlators are tree-level 
exact. We will relate this classicality to the fact that the sphere partition function, for any 
CFT, has a particularly simple c-dependence, and then discuss its bulk interpretation for 
holographic CFTs. Finally, we will discuss the generalizations of these statements at higher 
genus. 

The fact that correlators of descendants of the identity are independent of the theory, 
except its central charge, implies that we can compute them in any convenient theory with 
a variable central charge. One simple choice is the theory of c free bosons; by writing the 
relevant operators in terms of elementary fields, it is in principle straightforward to compute 
their correlators using free-field Wick contractions. This procedure is fairly tractable for 
calculating, for example, the four-point function of the stress tensor, but it rapidly becomes 
unwieldy when applied to higher-point functions or higher-level descendants, and for these 
calculations the methods described in the following subsections are far more efficient. 
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Nonetheless, the free-boson method gives a fast way to answer the important question of 
what powers of c appear in a given correlator. For example, the stress tensor is 


1 


T= -^ : aY'S.Y" : . 


(3.1) 


11=1 


An m-point function of stress tensors {T{zi) ■ ■ ■ T{zm)) includes indices pi,..., /im- The 2m 
X helds appearing can be contracted in various ways, linking the different T’s, and therefore 
the different p’s, to each other. For example, a contraction between and leads to 
a factor of ^^e connected part of the correlator, they are all linked in one group, 

so a non-zero contribution occurs only when all the p are equal: pi = ■ ■ ■ = pm- Hence the 
connected part of the correlator is linear in c, independent of m. Disconnected parts give 
higher powers of c; for example, the stress tensor four-point function has a term quadratic in 


c, from contractions in which the T are linked in two separate pairs (see (3.8) for the explicit 
form). 

All descendants of the identity can be written as normal-ordered products of derivatives 
of stress tensors. The connected part of any correlator—regardless of the number and type of 
operator—again comes from terms in which all the p’s are equal, and is therefore again linear 
in c. This also follows from the fact that the generating function for connected correlators 
of the stress tensor is the sphere free energy as a functional of the metric, which is simply 
c times the Liouville action. Thus, if we think of 1/c as playing the role of h, any CFT on 
the sphere is purely classical in this sector, in other words its correlators are effectively given 
entirely by tree-level contributions. To translate this into the usual held-theory language, if 
we normalize T by a factor of (so that its two-point function is 1), then a connected 
correlator with a total of P factors of T is proportional to just like a tree-level diagram 

with P external legs in a held theory with coupling constant 1/c. 

We now turn to holographic theories, where 1/c ~ Gjq is indeed the bulk coupling 
constant. An operator made of p stress tensors corresponds to a state containing p gravitons. 
This again leads to for a tree-level bulk process involving a total of P gravitons. 

From this point of view, the absence of loop corrections may seem mysterious, given that 
there certainly exist Witten diagrams in the bulk containing loops, which make non-zero 
contributions to such a correlator. However, in 3D gravity, all terms in the effective action 
which depend only on the metric (even those generated by loops of other helds) can be 
absorbed in the Einstein-Hilbert term |^. Hence, the full quantum ehective action for the 
metric is simply the classical action, with a renormalized value of the Newton constant. 
Since it is the renormalized Newton constant that enters in the relation 1/c = 2G7v/3i?AdS) 
when working in terms of c, the theory appears to be entirely classical]^ 


10 


This property is not directly related to the absence of propagating degrees of freedom in pure 3D gravity. 
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The arguments above, both in the held theory and in the bulk, depend crucially on 
the fact that we are working on the plane (or, more generally, on the sphere with any 


metric). On surfaces with non-zero genus, as we recalled in subsection 2.1, the free energy 


(or effective action) includes higher powers of 1/c, and these higher-order terms depend on 
the full operator content of the theory. So one would not expect correlators to be purely 
classical. Similarly, from a bulk point of view, gravitons and other particles can propagate 
in loops that wrap non-trivial cycles of the bulk, giving corrections that cannot be captured 
by a local effective action. 


Nonetheless, as noted below (2.15), at genus one, the vacuum free energy Fvac(T^) = 
— In Zvac(T^) does have the special property that it is one-loop exact, in other words contains 
only terms linear and constant in c (in any CFT). Zvac(T^) is dehned as the path integral 
with the insertion of the operator Pvac(^) that projects onto vacuum descendants on some 
fundamental cycle A of the torus. Derivatives of the free energy with respect to the metric 
give connected correlators of the form (Pvac(^)CliCl 2 ■ ■ ■)coin where the Oi are descendants 
of the identity. Such correlators are therefore also one-loop exact (contain only terms linear 


and constant in c). We will conhrm this property by explicit calculation in subsection 4.2 
below. 


3.2 Operator modes 

To begin, we will compute vacuum correlators of the form 


(O(«d)T(1)T(z)O(0)) 


(3.2) 


where {•) := {0| ■ |0). The operator O is allowed to be an arbitrary, non-holomorphic operator, 
not necessarily primary or quasi-primary. As is conventional, we leave its anti-holomorphic 
dependence implicit in what follows. We dehne mode expansions 


m = E 


^ 77.-|-2 


= E 


On 

^n-\-h 


(3.3) 


n G / 


n G ^ 


where the stress tensor modes obey the Virasoro algebra, 




mH-n,0 


(3.4) 


To demonstrate this, one can consider the correlator of four spin-s currents with s > 2. As we will show by 
example later in section 3.3.2| (see (3.41)-(3.42)), these correlators do not truncate in a 1/c expansion. This 
implies a non-trivial loop expansion for bulk four-point scattering of spin-s gauge fields in pure 3D higher 
spin gravity, even though these theories also lack propagating modes. 
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In terms of modes, the four-point function is 

{O{oo)Til)T{z)O{0)) = z-^ z-^{OhL-nLnO.H) . ( 3 . 5 ) 

n G Z 

To proceed, we break up the sum into the n = 0 mode term, and two sums over positive 
and negative integers (denoted and Z_, respectively). Using the fact that LoO-/i|0) = 
h(9_/i|0), the n = 0 mode contributes a term z~^h‘^J\fo, where Ao = {OhO-h) is the norm 
of O. Using the Virasoro algebra, and relabeling n —)■ —n, we can rewrite the sum over Z_ 
in terms of a sum over Z_|_ as 

z^(^{OhL_^KO.h) + (2nh + ^n{n^ - 1 )) ^^o) ( 3 . 6 ) 

n G Z_ n G Z+ 

= ^ z"{0hL.nLn0.h) + (2 ( 13 ^ + 2 ^( 13 ^) • 

The quantity in angle brackets is simply the squared norm ||L„C>_/i|0)|p. 

There is a further simplihcation of this sum: it truncates on account of vacuum invariance. 
Suppose (9 is a level N descendant of a primary held O' of holomorphic dimension H = h—N. 
Then O^h can be written as a linear combination of “lexicographically ordered” operators. 


L-m ■ ■ ■ L_n^.O_jj 


(3.7) 


where ni > n 2 > ■ ■ ■ > Uk and N = This implies that the sum in (3.6) truncates 

at n = A, because |0) = 0 for n > iV by dehnition of a primary. 

Taking this into account and adding (3.6) to the other pieces, the full correlator is 


(O(oo)T(l)T(z)O(0)) = 

/ 2 \ \ (3 8) 

^(z” •+ Z~''){OhL_nLnO-h) + (^2 (1 - z)4 ^^(1 _ zy ^ ' 

A pleasing feature of the expression in parenthesis is its manifest invariance under z ^ 1/z, 
which is simply invariance under crossing symmetry corresponding to exchange of the two 
stress tensors. 

If O is quasi-primary, then the n = 1 term of the sum vanishes. If O is primary, the 
entire sum vanishes. That leaves us with a very simple expression: 

C> primary: (C>(oo)T(l)T(^)C>(0)) = ATo . (3.9) 
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A useful check of (3.8) is to take O = T, which yields the stress tensor four-point function. 
Using the Virasoro algebra to compute (L 2 h_ 2 h 2 h_ 2 ) = c^/4, we hnd 


(T(oo)T(l)T(z)T(0)) = z-^{jfl + z^+ 


-I- 2c 


(l-z + z^) 


(3.10) 


(1 — z)‘^ j (1 — zy 

This agrees with previous results (e.g. (^). We may dehne the conformal cross-ratio as 

(3.11) 


^ 12^34 

X := - = ^ , 


^ 13^24 


in which case (3.10) has the correct form of a four-point function as determined by conformal 
symmetry: 

{T{oo)T{l)T{z)T{0)) = z-^T[x) . (3.12) 

We will reproduce this result in the next subsection in a more efficient wayp] 

One can easily generalize this analysis to correlators where T is replaced by a different 
operator. For simplicity, we consider the four-point function of two pairs of holomorphic 
quasi-primaries (9“ and (9^, of dimensions ha and respectively. (Their modes are dehned 


as in (3.3).) In this case, the resulting expression is 


{0“(oo)0‘(l)0‘(z)0“(0)) = 


(3.13) 




n=l 


n G 


To understand why the hrst sum truncates at n = we need to examine the OPE between 
quasi-primaries: in terms of modes. 


PI, 0^1 = Y. n-,K, ftj; K) + G“V+„,o 


m + ha — l 
hf) 1 


(3.14) 


is the Zamolodchikov metric, O'^ are also quasi-primary, 0“^^ are OPE coefficients, and 
the P(rn,n]ha,hb]hc) are known function^ encoding the contribution of the full global 
conformal family of All modes with n > —hf, annihilate the vacuum. This enables 
us to write |0) = [^ni^-h ]|0) ^ the OPE (3.14), combined with the 


unitarity bound h > 0, ensures that modes with n > ha give vanishing contribution. This 


^^Following our discussion in subsection 3.1 note that this could be derived for all c from a 3D gravity 
computation at large c, by thinking of T as a single-graviton state: the 0{P) part is the free-held Wick 
contraction, and the 0{c) part is the connected bulk correlator of four gravitons expressed in terms of the 
renormalized Newton constant. 

^^See e.g. equation (3.4) of 31 . All P{m,n-,ha,hi)-,hc) are hnite in unitary CFTs for operators of hnite 
dimension. 
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explains the upper bound in (3.13). 

Equation (3.13), while compact, is not particularly transparent. Even if is made of 
current modes alone, its modes may be given by inhnite sums over products of the L„, which 
are difficult to manipulate. More generally, the four-point function appears to depend on 
the full holomorphic operator content of the theory, due to the presence of the commuta¬ 
tor [0^,0^^]. In fact, this latter point belies the true structure of the result. We now 
demonstrate this explicitly as we turn to a much more powerful method of computation for 
correlators of holomorphic operators. 


3.3 The holomorphic bootstrap 

We will now describe a general method to compute the correlation functions of chiral oper¬ 
ators using crossing symmetry. We will see that any correlation function of chiral operators 
which obey a closed operator product algebra may be determined uniquely by a finite num¬ 
ber of three-point function coefficients. This is in contrast to the typical situation, where the 
OPE allows us to determine correlation functions only in terms of an inhnite sum over inter¬ 
mediate states. In many cases, such as for the correlation functions of Virasoro descendants 
of the identity, this leads to an extremely efficient computational algorithm. 

Let us recapitulate our conventions for chiral operators. We make no further reference to 
the mode notation of the previous subsection. We will consider a family of chiral operators 
Oa{z), with integer dimensions ha, and ha = 0. We will take the basis Oa to be quasi¬ 
primaries and assume that the Oa satisfy a closed OPE 


Oa{zi)Ob{z 2 ) ^ + (descendants) 

^19 


The two point functions 


and three-point functions 


{Oa{z,)Ob{z2)) 


Gab 


ha-\-hi) 

^12 


(a(o)a(^2)a(^3)) 




abc 


ha+h^-hc ha+hc-hi, h^+hc-ha 
^12 ^13 '^23 


are hxed, up to constants, by conformal invariance. 


(3.15) 


(3.16) 


(3.17) 
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3.3.1 Four-point functions: General case 

Conformal invariance constrains the four-point function to take the form 


{0aiz,)0tiz2)0,iz3)0diz,)) = 


ha-\-hij hc+hfi 


-12 


-34 


— I I — I I 

Zu 


Zl4 

Zl3 


(3.18) 


where we define the cross ratio x as in (3.11), 


X = 


2^12^34 

2 ^ 13^24 


1 — X = 


^14^23 

^13^24 


(3.19) 


We will use the notation H = ha, hat = ha — hi„ Zab = Za — Zb, etc. 


Our starting point is the observation that the four-point function (3.18) depends analyt 


ically on the Zi and has poles only when the points Zt coincide. Thus J^abcd is a meromorphic 
function of x with poles only at a; = 0,1, oo. So J^abcd is a rational function of x, which 
is uniquely completely determined (up to a constant piece) by its polar behaviour at these 
points. As we will see, this polar behaviour is fixed by only a finite number of three-point 
function coefficients. 

We begin by considering the expansion of J^abcd near a: —)■ 0. This can be found by 


inserting the OaOb and OcOd OPE into the four-point function (3.18). The result is a sum 


over intermediate operators Oe- The contributions from the descendant states of a given 
quasi-primary are given by a rigid (i.e. S'L(2,]R)) conformal block. The rigid conformal 
blocks were written in terms of hypergeometric functions in 35 . The result is 


J~ibcdi^^) ^ ^ {Cabe^cd ) ^ F{hg^ hab, hg -|- hcd, 2he, x) 


(3.20) 


From this we see that Fabcd is finite as x —)■ 0. The constant term as x —)■ 0 is given by the 
exchange of the identity operator, so 


FibcdijZ) GabGcd T . . . 


(3.21) 


where ... denotes terms that vanish as x —)■ 0. 

We now need to determine the polar behaviour near x —)■ 1 and x —>■ cx). To do this we 
will use the transformation properties of the four-point function under crossing symmetry. 
The crossing symmetry conditions can be derived by considering how the correlation function 


(3.18) transforms when the Zi are permuted. In particular, let us consider a permutation 
TT G 54 of four elements. We have 


{Ga(^Zi^Ob(^Z2')Oc(^Z^^Od(^Z4.^') (G7r(a) (■27r(l))G7r(()) (^7i-(2))G7r(c) (■^7r(3))f^7r(d) (^7r(4))) (3.22) 
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This relates J^abcd{x) to J^n{abcd){'^{x)), where the permutation vr acts on the cross-ratio as 


n{xj = — 

^7r(12)^7r(34) 


(3.23) 


One just needs to determine how the permutation vr acts on the prefactor in parenthesis in 


equation (3.18). Some permutations have 7r(x) = x; these give identities for the J^abcd{x) with 
fixed X. One can verify that these identities follow immediately from the conformal block 


expansion (3.20), using properties of the hypergeometric function identities and symmetries 


of the three-point function coefficients. Other permutations act on x, and give non-trivial 
information about four-point functions. In particular, the permutations tt = (14) and vr = 
(24) give the crossing symmetry equations 


J^abcd{x) = (-1)^X"“+"'^J-,,,,(1/X) 

= (-1)^X"^+"'*(1 - x)-^^-^^Fadcb{l - X) 


(3.24) 


These crossing equations strongly constrain the allowed form of the three-point function 
coefficients. Since J^abcdix) is finite as x —)■ 0, we see that J^abcd has a pole of order ha + 
at X —>■ cx) and a pole of order hb + he sX x ^ 1. 

We need to understand better the behaviour near these poles. To determine the behaviour 


near x —)■ oo we insert the conformal block expansion (3.20) into the first crossing symmetry 
equation to get 


J^abedix) = '^‘^F{he- hdb,K + hea-,2he-,l/x) (3.25) 

e 

ha^hd 

x" -f ... as X —)■ cxD . (3.26) 

n=l 

Here ... denotes terms that are finite at x —)■ cx). The important point is that, because the 
hypergeometric function is finite as x —>■ oo, the only terms that contribute to the pole are 
those with h^ < ha + In particular, the power series expansion of the hypergeometric 
function at x —)■ oo gives an explicit formula for the in terms of the three-point function 
coefficients CdbeC^a with he < ha + hd- We hnd 


a. 




ha+h^-n 

u ST- c 


C e 
ca 


he=0 


{he /iQ,+/id—/?-e—^“1“ hca)ha+h(i—he—r, 

{hd h(^ he Tl)\{2.he^ ha-^hd—he—n 


(3.27) 
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Similarly, near x —)■ 1 we have 


J^aU^) = {-If C'adeC'c6"(l - - Kd, K + Kb, 2^; 1 - 


hb+hc 


(3.28) 


/3n(l — x) ” + ... as X —)■ 1 


n=\ 


where ... denotes terms that are hnite as x —>■ 1. Again, the hypergeometric function has 
a simple power series expansion at x —)■ 1, giving an explicit formulas for the coefficients 
K in terms of the three-point function coefficients with he < K + he- The formula for the 
K is a bit more complicated than that for a„, since we must expand in powers of 

1 — X as well as the hypergeometric function, so we will not write it explicitly. However, the 
important point is that there is a completely explicit (albeit complicated) expression for the 
(3n in terms of the three-point function coefficients CadeCcF with he < K + K- 

The four-point function Fabcd is now completely fixed. It is the unique rational function 


of X which is finite everywhere except at 1 and oo, where its polar behaviour given by (3.26) 


and (3.28), and whose value at x = 0 is given by (3.21): 


ha+hd hb+hc 

•T ibcd{K) GabGcd T ^ ^ CXnX -|- ^ ^ (3n [(1 

n=l n=l 




(3.29) 


We see that Fabcd depends on a total oi H = ha + KF he + hd coefficients, /?„, which are 
determined by combinations of a finite number of three-point function coefficients. This is 
a consequence of crossing symmetry applied in a holomorphic setting; for non-holomorphic 
operators, there is no simple formula for a four-point function in terms of a finite number of 
operators. 

This has a remarkable consequence for the conformal bootstrap program, where crossing 
symmetry is used to place constraints on the three-point function coefficients. This is espe¬ 
cially true for chiral CFTs. In a typical CFT, the bootstrap results in equations involving an 
infinite number of three-point function coefficients, which can only be solved by truncating 
or approximating the crossing symmetry equations in some way. For a chiral CFT, the con¬ 
straints are all written in terms of a finite number of equations. For example, by comparing 


(3.20) with the expansion of (3.29) around x = 0 we can obtain explicit formulas for all of 
the coefficients CabeCcd^, for all e, in terms of the coefficients CdbeCca^ with he < ha + hd and 
CadeCcb^ with he < K + K. Of course, our results also apply to chiral operators in non-chiral 
CFTs. 


Moreover, we note that (3.29) is not the unique way of writing the four-point function. In 


writing (3.18) we chose to separate out a particular combination of Zij to define a meromor- 
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phic function. This choice led to a meromorphic function depending on H coefficients which 
were determined by three-point function coefficients CdbeCca^ with he < ha + hd and CadeCcb^ 
with he < hb + hc- Other ways of separating out a meromorphic function will lead to different 
expressions which in some cases may be more useful. For example one particular interesting 
way of imposing the crossing symmetry relations is to write the four-point function as 


Fabcdi,^^ 


ha~\~hf) — H/3 ha-\-hc—H /3 ha-\-hfi—H / 3 hff-\-hc — H/3 I3 hc-\-h^—H/3 


Z 


12 


2 : 


13 


2 : 


14 


23 


2:. 


24 


2:. 
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where 


J^abcdix) = ^"Fabcd{x) 


(3.30) 


(3.31) 


The function Fabcd is convenient because it treats the four points democratically - which 
makes the crossing symmetry equations very simple - but does so at the price of introducing 
a branch cut in Fabcd coming from the fractional powers of H/3. Fabcd has singularities of 
order H/3 at each of the three points a; = 0,1, cx); the crossing equations determine Fabcd to 
be 


lH/3\ 


Fabcdix) = _ a;)n-H/3) 


(3.32) 


n=0 


where the an,bn,Cn are determined by the three-point functions of operators with he < 
yH/3\. Note that, since n = 0,..., \_H/3\ we now have 3[Ff/3 -|- Ij coefficients to deter¬ 
mine. It is reasonably straightforward, though tedious, to write explicit expressions for these 
coefficients in terms of the three-point functions. The advantage of this approach is that 
it will, in principle, require the computation of fewer three-point function coefficients. For 
example, if the number of operators in the chiral algebra increases rapidly with dimension 
(as in the case of the Virasoro algebra) then this expansion would be much more efficient. 

3.3.2 Four-point functions: Identical operators 

Let us now simplify to the case where the four operators Oa are identical operators O of 
weight h, which is of interest for our computation of the higher genus partition function. 
In this case the above procedure simplifies considerably. The four-point function F{x) = 
Fabcdix) is a meromorphic function with poles only at a; = l,cx) which obeys the simplified 
crossing symmetry equation 


F{x) = x^^F{l/x) = x^'^i/i - xy^Fi/i - x) 


(3.33) 
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In fact, the space of such functions is a vector space of dimension 1 + [2h/3j. To see this, 
consider the function 


a{x) 


(1 — X + x'^y 
(1 — xy 


(3.34) 


which obeys (3.33) with h = 1. The function J^{x)a{x)~^ is invariant under the anharmonic 
group generated by x —>■ 1 — x and x —)■ 1/x. Moreover, this function is analytic everywhere 
on the Riemann sphere with the exception of a pole of order 2h at x = e”/^, along with a 
mirror image pole at x = These points are order-three hxed points of the anharmonic 

group, so when expanded around x = only cubic powers may appear. The function 


k{x) 


x^(i — xy 
(1 — X -I- x^)^ 


(3.35) 


is the unique meromorphic function invariant under the anharmonic group that has a pole 
of order 3 at x = We can therefore expand J^(x)a(x)“^ in integer powers of k, to 

obtain 


T(x) 


[2h/3\ 

n=0 


X 


2n 


(1 — X -I- X^) 


2h—3n 


(1 -x) 


2h—2n 


(3.36) 


To implement this way of computing four-point functions, we note that to determine the 
coefficients c„, we now simply expand this function in powers of x and use the OPE to 
determine these coefficients as products of three-point functions. 

It is instructive to phrase our conclusions in the language of modular functions. Equating 
X with the modular lambda function 


X = A(r) Ri — 128^ -|- , 


(3.37) 


where q = gives a map from Wlo ,4 (fhe moduli space of four marked points on the 

sphere) to (the moduli space of a torus). Accordingly, SL{2,Z) transformations of r 

induce anharmonic group transformations of x: specihcally, x —?• 1 — x and x —)■ 1/x are 
induced by the S and T ST transformations, respectively. The problem of Ending a function 
invariant under the anharmonic group therefore maps to finding a modular function, with 


desired polar structure in q determined by the poles in x via (3.37). Given the identification 


(3.37), our function fc(x) in (3.35) is just (256 times) the inverse of the J function: k{x) = 


256/J(r). So the construction of the four-point function J^(x) is literally identical to that 
of torus partition functions of holomorphic CFTs, as in [^. Likewise, (3.36) implies a 
Rademacher expansion for OPE coefficients of higher dimension operators. 
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We now treat some useful examples. For = 1 we have 


J^{x) = Co 


(1 — a; + 

(1 — xY 


(3.38) 


This is exactly the four-point function of a spin-1 current, j. The coefficient Cq = is 
determined by the hrst (trivial) OPE coefficient jjl, where k is the level of the current 
algebra. 

For h = 2 we get two possible functions, 




Co(l — X x'^)'^ -I- Ci(l — xYx‘^{l — X x"^) 

(1 — xY 


(3.39) 


This is the stress tensor four-point function. Matching the small x expansion with OPE 
coefficients of TTl and TTT, we hnd cq = c^/4 and ci = c(2 — c), where we used the 
canonical norm for the stress tensor, A/r = c/2 . This matches {TTTT) as computed in 

( ra . 

For h = 3 we have 




Co(l — X + xY^ + Cl(l ~ xYx‘^(l — X + x‘^Y + 0(1 — x) 




(1 


X) 


(3.40) 


This is the four-point function of a spin-3 current, call it hF, which was hrst worked out in 43 
in the context of CFTs with W 3 symmetry}^ Matching the small x expansion with OPE 
coefficients determines the Cj. In a theory with W 3 symmetry, the hrst three quasi-primary 
operators appearing in the exchange channel of (WWWW) are 1, T and the level-four quasi¬ 
primary A, which is the normal-ordered product of T with a derivative term subtracted: 


A := (TT)- d'^T . 

^ ^ 10 


(3.41) 


This operator has norm TVa = c(5c -f 22)/10. The OPE coefficient WWW vanishes, as it 
does for W being any odd-spin chiral quasi-primary. Computing OPE coefficients using the 


Virasoro algebra and matching to the small-a; expansion of (3.40), we hnd 


Co — , Cl — A/^ 


w 


6(3 


C2 


= M 


,2 I 3(5c 2 - 71c - 102) 
w 


c(5c + 22) 


(3.42) 


where {W 3 W- 3 ) = Afw oc c is the norm of W^ As explained in section 3.1, {WWWW) 


believe our result actually corrects a sign error in |43| : the parameter /r there should be a sum, not 
a difference, of two terms. 

the chiral algebra contains a spin-4 current too, as in the case of the Woo [A] algebra appearing in 
the context of higher spin AdS/CFT [^, this current will also appear at level four with nonzero OPE 












does not truncate in a 1 /c expansion. 
For = 4 we get 


T{x)^ 


Co(l — X + + Ci(l — x^x'^il — X + x'^Y + ^ 2(1 — xYx^il — X + x'^Y 


(3.43) 


(1 — xY 

An example of an = 4 chiral operator whose four-point function we will need in the sewing 


construction of Zy^c is the quasi-primary A introduced in (3.41). Again computing the q by 
matching to OPE coefficients, we hnd 


c, = Ar'i, 


,2 , 4(125c2 + 590c + 3704)\ 


5c(5c + 22) 




(3.44) 


In closing, one interesting comment is that certain general consequences can be immedi¬ 


ately read off from (3.36). For example, when h is not a multiple of 3 every contribution to 


the four-point function J^{x) includes a factor oi 1 — x + x'^. Thus if h is not a multiple of 
3, the four-point function vanishes when x = ±e'^*/^. 


4 Genus-two partition functions 


We are now ready to compute Zvac at genus two, which captures the contribution of the 
Virasoro vacuum module to the partition function of a CFT on a genus-two Riemann surface. 
We begin this section by reviewing the Schottky uniformization of generic genus-^f Riemann 
surfaces and the sewing construction of the genus-gf partition function. We then turn to 
the actual computation of Zvac at genus two using sphere four-point functions of low-lying 
Virasoro vacuum descendants. The hnal result can be found by substituting the results of 
subsection 4.2 into equation (4.12). We focus on the holomorphic part of Zvac henceforth. 


4.1 Schottky uniformization and the partition function 


A non-singular genus-^f Riemann surface can be constructed by cutting out 2g disks on the 
Riemann sphere and identifying pairs of boundary circles to form g handles. The Schottky 
uniformization of the Riemann surface entails the identihcation of pairs of circles through 
Mobius transformations 7 *, i G {1,2, •• • ,g}, which are elements of PSL{2,C). The maps 
7 j form the generators of the Schottky group, F. There are three parameters {ai,ri,pi} 
associated with the handle: {ai,ri) are the locations of the centers of the boundary 
circles, and pi determine the width of the handles. A global conformal transformation can 


coefficient, and will change the value of C 2 . 
algebra; one instead finds C 2 = 3A/'^(A^(25c^ 


This agrees with a previous result 


45 


This generalization is simple to compute using the Woo [A] 
- 115c + 546) - lOOc^ - 740c - 7464)/(5c(5c -f 22)(A^ - 4)). 
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fix the positions of three boundary circles on the sphere and thus a genus-^f Riemann surface 
has 3g — 3 complex moduli. 

We consider the Schottky uniformization of a genus-gf Riemann surface following the 
conventions of (see their Appendix C), where the locations of each pair of identihed 
circles are given by the Mobius transformation 




TjZ + aj 
^ + 1 ’ 


(4.1) 


where 7ai,ri(0) = Oj and ■jai,ri{oo) = r^. The generators of the Schottky group are given in 
terms of this map as 

li = la^,nlpna!,n^ (4-2) 


where 'jpX^) = Pi^- We note that identihed circles have opposite orientations: for the 
pair the two boundary circles are given by the maps Ci = '^ai,rilRiC and C-i = , 

where C is the unit circle at the origin, Ri and R-i are the radii of Ci and C-i respectively, 
and 7 is the inverse map 


7(z) = 


(4.3) 


The product of the radii of the two circles is Ri R-i = Pi. We refer the reader to Appendix 
C of 40 for more details on the Schottky parametrization. 

The partition function of a genus-gf Riemann surface uniformised by the Schottky group 


is given by the following power series expansion in pi 40 


15 


^9=^1 Pi" ---Pg" ^hu...,h,{a3,...,ag,r2,...,rg). (4.4) 

hi,-- , hg 


In the sewing construction, a handle is replaced by the boundary states inserted at the 
centers of the two disks. The functions Chi,--,hg are 25 f-point functions on the Riemann 
sphere and R is the conformal dimension of the operators inserted at the pair of disks. 
A schematic picture of the sewing construction is shown in hgure In the above equation 
we have hxed the positions ai = 0, ri = oo, and 02 = 1. The 25 f-point functions Ch^^...^hg-i 
whose ingredients we will explain in the next paragraphs, are sums over products of vertex 
operators of the form 


Chi,-- ,hg — 




E n CA ( n r,) V"(*, a,) 


2=1 


(4.5) 


Actually, the formula (4.41 just gives the partition function up to a factor of the form e in other 


words, in an expansion of the free energy F := —InZg in 1/c, it only gives the order C and higher terms. 
The order-c^ term cFq depends on the full metric on the Riemann surface, not just the complex structure. 
Its calculation within the context of the sewing construction is explained in appendix 
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where G is the Zamolodchikov metric dehned below in (4.20), and "H/i- is the Hilbert space 
of states of dimension hi. 

These expressions for the vertex operators should be understood as follows. Under any 


Mobius transformation 7 ( 2 ;), the vertex operator V{(j){z)) transforms as 46 

u|^t/(7(z))0,7(z)^ = 


where 


'y'iz) = ^ 
dz^ 


/// N (Pi 


For the Mobius transformation (4.1) we have 


y (z) = ^^ 7 " (z) = -2 — _ — 


' . [ai-Ti) 

la,,nl G) = 




7a„r,7) ( 2 ) = -2 


(a* - U) 

(;2 + 1)3' 


(4.6) 


(4.7) 


(4.8) 


The “in” and “out” vertex operators in the sewing construction (4.5) then transform as 


ryy,a,) = v(^[/(ja.,rAz = 0))y,7a„.,(;^ = 0)^ = (n - a,)^°e-^V,(a,), (4.9) 


and 

= P))Phlai,nl{z = 0 )^ 


(_l)^o(^._^.)iog 'f^iy.(rj). (4.10) 


For i = 1, he. for the handle with the two ends at (0, oo), one can perform an extra Mobius 
transformation under which the two maps at zero and inhnity become the identity and the 
inverse map, respectively. This is described in the next subsection. We note that if and 
xjji are quasi-primaries, then the vertex operators are given by 


V^"'{(j)qp,ai) = (fi - ai)^-(’<ip(j)qp{ai), 
V°^\'il)qp,ri) = (-l)^’^'3P(ri - ai)^'i’i^'tpqp{ri). 


4.1.1 Genus two 

We now specialize to genus-two Riemann surfaces. The partition function is given by 

(4.12) 


OO 

hi,h2=0 
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where we have defined r 2 = x. Using (4.5), the functions Chi,h 2 {^) are found to be 


These two formulae apply in general. For our purposes of computing Zvac, we only allow 
Virasoro descendants of the identity to be inserted at the boundary circles of the handles 
as in figure]^ Henceforth, we refer to (4.12) with the understanding that we compute Zvac 
specifically. 


Let us now define the vertex operators needed in (4.13), starting with those at (0, oo). 


The functions Chj^^h 2 {x) ai'e invariant under the map 7a,,ri —t 'Jai,ri1t, where 'yt(z) = tz, 
t E C*. For the i = 1 handle with its two ends located at ai = 0 and ri = oo, we consider a 
Mobius transformation of the form 7ai,ri7i/ri and find 


7ai,ri7^ - ^ , A 

ri 


= Z, 


ri—>oo 


7ai,ri7^7 = —X 
’"1 z + — 


ri—>-oo Z 


(4.14) 


This therefore gives the identity map for ai = 0 and the inverse map for ri = cx). The vertex 
operator at the origin is simply 


V”(.).1,0) = V'(.^1,0) = *(0). 


(4.15) 


The vertex operator at infinity follows from using 


1 




(4.16) 


which yields 


oo) = oo^ = Im (- 1 )^°z 2 Lo (4.17) 


For the handle with vertices at ( 02 ,^’ 2 ) = (l,x), we use (4.9)-(4.10) to read off 


r"(02,l) = (x-l)^°e-^V2(l), 


(4.18) 

(4.19) 
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The Zamolodchikov metric is defined in terms of the in and out vertex operators a: 


G« = ^V“‘(^i,oo)r”(.#..0)y 


(4.20) 


We choose an orthogonal basis of states at each level by diagonalizing the Gram matrix. The 
Zamolodchikov metric is thus diagonal and the ingoing and outgoing vertex operators are 
the same up to Mobius transformations. Consequently, we can (and will) define the norm of 
the states as 

To summarize, the following is the prescription for constructing the genus-two partition 


function: insert vertex operators (4.15), (4.17)-(4.19) and the Zamolodchikov metric (4.20) 


into (4.13) to evaluate Ch-^^h 2 i^): and sum over these using (4.12) 


4.2 Four-point functions Ch^^h 2 {x) 


We will momentarily compute some of the functions Ch^^h 2 {.x) defined in (4.13). Before doing 
so, it is useful to elucidate some of their general properties. 

First, these functions are symmetric under the exchange of the positions of the two 
handles: 

= Ch2M{^) (4-21) 


When hi = 0 or h 2 = 0, 


Co,h{x) = Ch,o{x) = d{h) , 


(4.22) 


where d{h) is the degeneracy of operators at level h. This follows from the definition of the 
vertex operators and of the Chi^h 2 {x) themselves. It can also be understood intuitively: re¬ 
placing either handle with two insertions of the identity reduces the genus-two partition func¬ 
tion to the torus partition function, the holomorphic half of which is Trp^“ = d(h)p^|^ 


Thus, (4.12) implies (4.22). 


Now we consider the x-dependence of Chi,h 2 {x). For general x, they obey 


Ghi,h2{x) (^/^) 


(4.23) 


which is a consequence of modularity with respect to Sp{4,7j). Taking various limits in x 
corresponds to taking OPE limits of the (dressed) four-point functions defining Chi,h 2 (x). 


^®We note that our formulae (4.10) and (4.17) contain an extra factor of (—1)'^° comparing to the formulae 
in Appendix C of 40 . The reason is that we choose a different convention than that of 40 . In our convention 


G is the Zamolodchikov metric whereas in j40j their metric G is a metric on the space of states which is 
related to G via G^^ = G(_i)i.o 0 ,^. 

^^Note the c-independence of this quantity. When summing over the vacuum module only, the dimensions 
of each state are fixed by conformal symmetry, and hence unrenormalized. This is the CFT statement of the 
one-loop exactness of the pure gravity partition function on a solid torus. 
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The simplest one is x —)■ 1, which describes the fusion of two ends of the same handle. In 
this case, 

= d{h) X dih) . (4.24) 


This again follows by dehnition, and is necessary for the partition function (4.12) to factorize 
in the separating degeneration limit. 

More subtle are the equivalent OPE limits x —?• 0 and x —>■ oo, which describes the fusion 
of two ends of different handles. These limits are singular. In Appendix]^ we show that in 
a 1/c expansion, the leading powers of x —)■ 0 that appear are correlated with powers of 1/c 
as follows: 


lim \imCh,,h-hiix) ~ 0(x '*) + -0(x 4 1 

£—>•00 x^O C \ ' 

\n=2 


(4.25) 


We have dehned h 2 = h — hi, and are assuming h > hi > 0 because C'o,/i(x) is constant. We 
are ignoring hi- and h-dependent coefficients at each order, and displaying only the leading 
singular behavior at each order in 1/c. The last term means that at (P(l/c^) and all orders 
beyond, the leading divergence scales as 0(x“^’''^). 

We now proceed to compute Chi,h 2 {x) explicitly for low values of hi and h 2 . A word on 
notation: henceforth, we denote the set of operators at level h above the ground state as 
where i = 1,2 ,..., d{h). 


4.2.1 

In this case, the identity operator propagates through one of the handles, so the four-point 
functions reduce to two-point functions 




O.,o(x)= ^ G4(l/“‘(-iA,oQ)r"(.#.,0)), 




(4.26) 


As discussed earlier, Ch^^x) = d{h)-, from the dehnition of in (4.20), this is obviously 


true. It is less obvious that Co,/i(x) = d{h) by looking at (4.26). Hence, we hnd it instructive 
to compute two CQ^h{x) in detail, to illustrate how to use the method of outlined above. 
In the hrst example we consider a quasi-primary state and in the second example we con¬ 
sider a secondary state. The latter is particularly useful, as secondary operators transform 
nontrivially under the Mobius transformations, so care must be taken in computing their 
correlation functions. 
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At level 2 {h = 2) there is only one state, the stress tensor T = L_ 2 | 0 ), which is a 
quasi-primary with norm A/r = From (4.18)-(4.19), the vertex operators at x and 1 are 


y«“*(T,x) = (x- 1)2 T(x), 
1/*"(T,1) = (x- 1)2T(1). 


Then we have, as expected. 


CoA^) = ( V°^\T,x) V^^{T, 1 ) > = 1 


(4.27) 


(4.28) 


where we used the stress tensor two-point function 


{T{x)T(l)) 


c 1 
2 (x - 1)4 ■ 


(4.29) 


At level 3, there is again one state O 3 = dT = L_ 3 | 0 ), this time a secondary state with 
norm Aog = 2c. The vertex operators at x and 1 are now, from (4.18)-(4.19), 


1 /°“'(C> 3 , x) = -(x - lfdT{x) - 4(x - l)^T(x), 
1 /*^(C> 3 ,1) = (x - 1)^5T(1) - 4(x - 1)2 t(1). 


(4.30) 


Using (4.29) we again have, as expected. 


C'o,3(a;) = MA,(v^\0^,x) ^"( 03 , 1 )^ 

= ^^ ( - (x - 1)^ dT{x) - 4(x - 1)2 T(x)) (^(x - l)^aT(l) - 4 (x - 1)2T(1)) ^ 

= 1. (4.31) 


It is also useful to see the vertex operator at inhnity, which carries non-trivial dressing on 
account of O^, being a secondary operator: 

U°“‘(C> 3 ,cx)) = lim uf(-l)^°^ 2 Log-^Li( 9 g^^A = {-z^dT{z)-Az^T{z)\. (4.32) 

We next move on to the computation of Ch^^h 2 A) 7 ^ 0 and ^2 7 ^ 0. We compute the 

requisite four-point functions using the methods described in section The transformation 
properties of the vertex operators are evaluated following the same procedure shown in the 
above examples; accordingly, the presentation here is streamlined, with some details of the 
vertex operator transformations relegated to Appendix There, we also list the operators 
and their norms through level six of the vacuum module. 
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4.2.2 


First, consider the four-point function with ^2 = 2 which corresponds to the four-point 
function of four stress tensors, each dressed with appropriate factors of x and 2 ; (recall 


equation 4.13). We hnd that 

1 


lim { j-TW (x-lYT(x) (x-iyT{l) T{0) 

(I) 


= 1 + (i - 1)* + 


4 ^ (x — 1 )"^^ ^ 8 (x — 1 )^ (1 — X + x^) 


x^ 


H— 
c 


x^ 


(4.33) 


We note that C 2 , 2 (x) is manifestly symmetric under x 1/x, as required. 

We next compute 6 * 3 , 2 (a;)- This can be done by taking the derivative of 6 * 2 , 2 (a^) with 
respect to the sphere coordinates at the insertion points of O 3 , or by direct computation of 


the four-point function using the Virasoro mode expansion formula (3.8). We obtain 

6 * 3,2 (a^) = ^ (4 — 16x -I- 24x^ — 16a:^ -f ix'^ + x^ + 4x® — 16x'^ + 24x® — 16a:® -|- 4a;^® j 
1 2(x- 1)2 


x-^ 


4 -I- X — 4x -f X -f 4x ), 


(4.34) 


We used (4.27) and (4.32) to write down the necessary vertex operators. One can explicitly 


check that 6 * 3,2 (x) = 6 * 2,3 (x), as required. 

We next move to level four and evaluate 6 * 4,2 (x). There are two orthogonal states at level 
four, 

(4.35) 


=A={L 


.4 |0), of = L_4|0) , 


where A is the quasi-primary hrst recalled in (3.41), and Of = d'^T/2 is secondary with 


norm J\i.^( 2 ) 


,( 2 ) — 5c. We obtain 
2 


6 * 4,2 (a:) ^ ~ 20a: + dlx^ — 24x^ -|- llx"^ — 4x^ -|- 3x® 

— 4x^ -I- llx® — 24x® -f 31x^° — 20x^^ -|- 5 x ^2 
1 4 (x — 1)2 


(4.36) 


4 — 3x + lOx — 14x -f lOx — 3x -f 4x 


c X* 

Note that the hnite expansion in 1/c is at hrst surprising because A/yC^ has an inhnite 1/c 
expansion; there are non-trivial cancellations between the two correlators. We will explain 
this in a moment. 

For the remaining computations, we will be briefer. (We remind the reader of Appendix 
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[^containing further details.) At level five, we find 


Cg 2 ( 0 ;) = 4 f20 - 80x + l2Ax^ - 95x3 ^ 

’ x‘ V 

+4x^ - lOx® + 40x^° - 95x3^ + 124x^2 - 80x^3 + 20x3^ 
1 4(x - If 


(4.37) 


c x° 

At level six, we find 


6 — 6x + 7x + 5x — 14x + 5x + 7x — 6x + 6x 


Cg 2 (x) = ^ (^35 - 140x + 220x2 - 172x3 + 67x^ - Sx^ + 2x® - 8x^ + 12x3 (4.38) 

-8x^ + 2x30 - 8x33 ^ g7^i2 _ ^72x33 + 220x3^ - 140x3*3 + 35 x3o) 

+ - (s — 7x + 10x2 — 8x3 _|_ _|_ ^g^6 _ g^7 

C X° \ 

d-lOx® — 7x0 _|_ 3^10j _ 


We observe that all the functions Ch^ 2 (x) we have computed so far contain a term pro¬ 
portional to 1/c and a term constant in c. In fact, this is true for all h. The reason is that 
the sum over h of Ch, 2 {x) is related to the two-point function of the stress tensor on the 
torus with the insertion of a vacuum projector: 

(P„or(l)T(x)> = (4.39) 


(where pi = This two-point function can in turn be obtained by differentiating the 

vacuum free energy with respect to the metric. This free energy was shown at 

the end of subsection 3T to be one-loop exact, in other words to contain only terms linear 
and constant in c. Hence, for all h, Ch, 2 {x) must contain only terms constant and inversely 
proportional to c. 


4.2.3 

We next compute the functions Ch^^hfx) with hi > 3 and /i 2 = 3. The hrst function is 
Cs^six) which corresponds to the four-point function of four cIT’s. This correlation function 
can be evaluated by taking the derivatives of C 2 , 2 {x) with respect to the sphere coordinates 
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at the location of the four operators, or by direct computation. We hud 


C3 3(0;) = L (25- IlOx + - 164 a;^ + 71x^ - Ux^ ( 4 . 40 ) 

’ x° V 

+3a;® - 14a:^ + 71x® - 164a:® + 191a:^° - 

+ - ——f 18 — 6 x + — 8 a:^ + x'^ — Qx^ + 18a:®), 

c x^ \ / 

Next, for ( 74 , 3 ( 0 :) we hud 

Ci 3 ( 0 ;) = 4 ("45 - 210a: + SOga:^ - 3960 :^ + 219a:^ - 66 a:® + 9a:® (4.41) 

’ x‘ V 

W + 9a:® - 66 a:® + 219a:^® - 396a:^^ + 399a:^® - 210x^® + 45a:^^) 

+ - ~ (qA- 83a: + 68 a:® + 3a:® - 56a:^ + 3a:® + 68 a:® - 83x^ + 64a:®), 

c x^ \ J 

and for (^ 5 , 3 ( 0 :) we hnd 

Jx) = 4 (245 - 1190a: + 2380a:® - 2526a:® + 1530x^ (4.42) 

’ x° V 

-530a:® + 100a:® - lOa:^ + 4a:® - lOx® + 100a:®® 

-530a:®® + 1530x®® - 2526x®® + 2380a:®^ - 1190a:®® + 245a:®®) 

1 (X — 1)^ / 

+ - ^^ (150 - 264a: + 201a:® - 32x® + 3x^ - 56a:® + 3a:® 
c x° \ 

-32a:^ + 201a:® - 264a:® + 150x®®). 


Again, all of the functions Ch, 3 {x) evaluated so far truncate at order 1/c in a large-c 
expansion. This is because all four-point functions in Ch, 3 {x) can be computed by taking 
derivatives of stress tensors in Ch, 2 {x), which do not affect the c-dependence of the correlators. 


4.2.4 C^,4(x) 

The function (^ 4 , 4 ( 0 :) is a linear combination of four four-point functions: 


* 5=1 


CiA(x) = 5^ ( r”“(04,, CO) V“‘(04j, I) V"(04,4, 1) V"(04,., 0) ) . (4.43) 


By the same argument below (4.37) and (4.42), the terms which contain at least one pair of 


the secondary operator O^®^ = d‘^T/2 truncate at order 1/c in a large-c expansion. The only 
term in ( 74 , 4 ( 0 :) which could potentially contribute at higher orders in 1 /c is the four-point 


function of four quasi-primaries A, dehned in (3.41). Let us focus on this contribution, which 
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we call C 4 ^ 4 |A(a;). 

Using the definitions of the vertex operators simply yields 

— 1)® lim z®/A( 2 ;) A(a:) A(l) A(0)\ . 

2:^00 \ / 


(4.44) 


The norm of A was given in (B.4). Substituting this and the results (3.43)-(3.44) for the 


four-point function obtained via the holomorphic bootstrap, we find 

(x — 1)^(1 — X -I- 


„ I , . (1 - X -t- x^)® /32 

C'4,4|a(x) = ^^^ + ( — 


X° \ C / X' 

,2^ 1 ^4 


6 


2\2 


4 (3704 + 590c + 125c2) (x - 1)^(1 - x + x^) 


5c (22 + 5c) 


x^ 


(4.45) 


Crucially to what follows, we observe that C 4 ^ 4 |a(x) contributes to an infinite expansion in 


1/c. This comes entirely from the inverse norms in (4.44) 


The collection of C/jj^/i 2 (x) computed in this subsection, plugged into the partition func¬ 


tion (4.12), forms one of our main computational results: namely, the first several terms in 


the Virasoro vacuum module contribution to the partition function of an arbitrary CFT on 
a genus-two Riemann surface, in the regime of Schottky parameters Pi,P 2 ^ 1- 


5 Free energy at large central charge, 3D gravity, and 
Renyi entropies 


Having derived the first handful of terms in equation (4.12) for all c, it is trivial to expand it at 


large c. As we have discussed, this large c expansion may be interpreted as the semiclassical 
expansion of the pure 3D quantum gravity partition function around genus-two handlebody 
geometries with conformal boundary S specified by Schottky parameters {pi,p 2 ,x}. In 


subsection |5.1[ we present some of our main results. First, we provide explicit contributions of 
the Virasoro vacuum representation to the CFT free energy at all orders in 1/c, corresponding 
to all-loop free energies in the gravitational loop expansion. We also show that at least in 
the perturbative regime pi,p 2 1, the loop expansion does not truncate except when S is 
the union of two tori. 


We then proceed to subsections 5.2 and 5.3, where we expand our general result (4.12) 


near two symmetric points in the genus-two moduli space: the replica surface ^ 2,3 used to 
compute the Renyi entropy S 3 for two disjoint intervals in vacuum, and the point corre¬ 


sponding to the separating degeneration limit x = 1. Our results will extend those of 33 
and (28], respectively. 
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5.1 All-loop results in 3D quantum gravity 

We consider the 1/c expansion of the vacunm free energy, Fyac = — logZvac: 

OO 

^vac = C^“^Fvac; i (5.1) 

i=0 


where i denotes the loop order. Its moduli-dependence is kept implicit. We can read off the 
loop corrections Fvac;£ from (4.12) npon expanding the Chi,h 2 {^) 1/c. We likewise expand 

these as 

OO 

. (5.2) 

£=l 

To begin, note that in the small {pi,P 2 ) expansion in which we work, both the one- and 
two-loop free energies are nonzero. This follows from the explicit resnlts in section and 
from the Renyi entropy compntation (2.15), but also from our general exposition of c-scaling 


of identity module correlators in section 3.1 


More interesting is the question of whether there are higher-loop terms. For i > 2, the 
Chi,h 2 {x) that we have computed all obey Chi^h 2 -,£{x) = 0 except for computed in 

(4.45), which clearly has an inhnite expansion. Accordingly, the leading contribution to the 
three-loop free energy Fvac ;3 in a small (pi,P 2 ) expansion can, and does, appear at O^pfp^)- 


Fva.c-,3 — 


ptpt ^ {C2,2-,2{x)f^ + 0 {p\pI) . 


There is no cancellation: instead, our results yield 




vac; 3 


4 4 13312(x 
P 1 P 2 - 


iy{i-x + xy 

2bx^ 


+ 0{ptpl) . 


(5.3) 


(5.4) 


For a; G M, this only vanishes at a; = 1. This point in moduli space corresponds to the strict 
separating degeneration limit. As the torus free energy is known to truncate at 0{c^) in a 
large c expansion [^, the fact that FVac ;3 = 0 when x = 1 is required by consistency. The 
interesting result proven here is that Avac ;3 is nonzero everywhere else on the real linep^ 

In fact, the inhnite 1/c expansion of < 74 , 4 (x) and the hnite 1/c expansion of C 2 , 2 ix) 
together imply that the f'-loop free energy Avac;£ is nonzero for all i, at least in a small 
{pi,P 2 ) expansion. The reason is simply that the term at 0 (PiP 2 ) cannot be cancelled by 
higher order terms in pi and p 2 . For i > 3—that is, at 0{l/c^) and beyond—FVac;r is given 


^®Note that (5.3) and (5.4) also hold when pi = p 2 , because the only contribution at C>(l/c^) to Chi,h 2 {x) 
for hi + h 2 = 8 comes from ( 74 ^ 4 (x). On the other hand, if x is a function of pi and p 2 , higher order terms in 
(5.4) are not necessarily suppressed. We will encounter such a situation in our discussion of Renyi entropy. 
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to leading order in pi , p 2 as 


F„e>3 = PiPi <^4,4; 03 ( 2 :) + 0(p;p|) 


(5.5) 


with 


G 


4,4; £>3 


(x) = 


X 


1)^(1-x + x2)2 /4(3704 + 590c+125G) 


x^ 


5c(5c + 22) 


(5.6) 


determined by (4.45). From a technical standpoint, this non-trivial 1/c expansion arises 
from the inverse norms appearing in the sewing construction, cancelling the norms in the 
correlator (AAAA). 


Interpreted as a CFT result, (5.5) is an exact expression for the contribution of the 


Virasoro vacuum module to the genus-two free energy of any family of CFTs that admits 


a 1/c expansion. Interpreted as a pure gravity result, (5.5) is an explicit formula for all¬ 


loop free energies on genus-two handlebodies. The loop counting parameter in the bulk is 
Gn = 3i?Ads/2c. In contrast to the one-loop exactness at genus one, the genus-two free 
energy is not exact at any loop order. 

Strictly speaking, we have so far established that the semiclassical expansion does not 


truncate for any real x 7 ^ 1. What about complex x? In particular, the all-Ioop terms (5.4) 
and (5.6) clearly vanish at x = the complex roots of 1 — x - 1 - x^ = 0. This follows 


from the same property of the four-point function of A, as discussed below equation (3.44). 


But this is a special feature of correlators of identical operators with h/3 ^ Z, so it will 
not persist to higher orders in the sewing expansion. For instance, G 4 ^ 5 ;^> 2 (e^®®®/^) 7 ^ 0, and 
likewise at all higher levels. Therefore, we have shown the following statement: perturbatively 
in {pi,P 2 ), the loop expansion does not truncate for any genus-two handlebodies except at the 
separating degeneration point. 

We note that at hxed order in pi and p 2 , the 1/c expansion for Zvac converges. This 
does not necessarily imply, however, that at a hxed point in moduli space (i.e. for hxed 
values of pi and P 2 ) the 1/c expansion converges. Indeed, since there are presumably other 
saddle point contributions to the path integral (coming from bulk handlebodies with diherent 
topology) one might expect that the series expansion of is asymptotic in 1/c. However, 
in the genus one case the series expansion for Z^ac converges—in fact, it truncates at order 
c°. Thus it is an interesting open question whether the 1/c expansion for Z^ac converges at 
higher genus. 

Finally, let us comment on the holographic interpretation of our result for the one-loop 
partition function, Zvac;i- This can be viewed as a computation of the holomorphic half of 
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the graviton handlebody determinant d® 


Zf- = J] J] 


7 G7^ n=2 




(5.7) 


The novelty of our computation is that we work in the regime of pi,p 2 ^ 1 hut for arbitrary 


X. This regime has not been probed directly in existing computations of (5.7). In 15 , (5.7) 


was computed for handlebodies asymptotic to replica manifold for two-interval Renyi entropy 
in a short interval expansion; this has only a single modulus and requires Pi = P 2 ^ 1 and 


a; ^ 1. In 28 , (5.7) was computed near the separating degeneration limit where S becomes 
the union of two tori, which requires x ~ 1 for arbitrary (pi,P 2 )- 


5.1.1 Higher spin theories 

So far, we have restricted to the pure Virasoro sector of the CFT. The meaning and calcula¬ 
tion of Zvac are conceptually unmodihed in the presence of higher spin currents. Along with 
the stress tensor, these Virasoro primaries live in the vacuum representation of an extended 
conformal symmetry, typically a W algebra. In the computation of Zvac by sewing, we now 
allow these currents and their normal ordered products to propagate through the handles. 
The resulting Zvac is again of the form (4.12), only with different coefficients Ch-^,h 2 {.^)- 
The holographic dual of Zvac in the presence of higher spin symmetry is the perturbative 
partition function of pure 3D higher spin gravity. A bulk Chern-Simons theory with connec¬ 
tions valued in two copies of a Lie algebra Q describes the vacuum sector of a CFT whose W 
algebra is the Drinfeld-Sokolov reduction of ^ 47 . Accordingly, the 1/c expansion of Zvac 


for such a CFT yields the semiclassical loop expansion of the Q x Q Chern-Simons higher 
spin theory. 

As a simple example, consider a CFT with symmetry, which contains a single higher 
spin current of spin three, W. Its presence will modify most of the Ch^^h^{x) coefficients, 
starting with Cz, 2 {x) and C 3 , 3 (x)j^ We can easily compute these using the correlators of 
section]^ The interesting term is (73,3(x). Denoting the contribution to (73,3(x) from the W 


19n 


The first product in (5.7) runs over primitive elements 7 G F C T, defined as those elements that cannot 


be written as 7 = /3"* for /3 G T and m > 1. The eigenvalues of 7 are eig( 7 ) = 97 and we do not count 7 
and 7 “^ as distinct elements. 

^°The generating function of quasi-primaries containing at least one W current is given in Appendix B 


of 48 . 
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current four-point function as 6wC^^^{x)^ we find, using (3.40) and (3.42), 


5wC^,z{x) := Mw^{x - 1 )® lim z\W{z)W{x)W{l)W{Q)) 

Z^OO 

(1 — a: -I- 6(3 — c) (1 — a:)^(l — a; -|- a;^)^ 

x^ c x'^ 

3(5c 2 - 71c- 102) (1 -a;)^ 

^ c(5c + 22) ■ 


(5.8) 


Expanded at large c, this yields an infinite series of loop corrections to the free energy of 
Q = SL{3) higher spin gravity: 


+ 0{plpt). (5.9) 

This is nonzero for all x 7 ^ 1, so we conclude that the loop expansion does not truncate 
away from the separating degeneration limit for small (pi,P 2 )- This is true for all higher spin 
algebras Q. We will return to the topic of higher spin theories in the Discussion. 


IpSL(3) _ 3 3 (1 ~ ’ 

Ka,c -,£>2 — P1P2 ^2 


3(5c 2 - 71c - 102) 
c(5c + 22) 


5.2 Renyi entropies 


As discussed in section there are three Renyi entropies that involve genus-two replica 
manifolds without punctures (i.e. for CFTs not in excited states). These are the N = 2,n = 3 
and N = 3,n = 2 Renyi entropies for a CFT on the plane, and = l,n = 2 for a CFT 
on the torus. We mostly focus on the N = 2 case, with replica manifold ^ 2 , 3 - Our results 
in section 5.1 are sufficient to rule out the truncation of the 1/c expansion of TVac even in 


'•2,3 


the case of the replica manifold S = 
detail; the hnal results can be found in (5.17) and (5.18). 


introduced in section 2.1 We now exhibit this in 


5.2.1 Two intervals on the plane 

Our goal is to express the free energy in terms of the coordinate y parameterizing the interval 
spacing, dehned in section To do so, we need only to express the Schottky coordinates 
{pi,P 2 , x} in terms of y. One way to proceed is by using the period matrix D(|/) for the replica 
manifold .^ 2 , 3 , which is known 38 . Thus, we will perform the map {pi,p 2 ,x} h-)■ {(?p( 2 /)}. 


where qijiy) = exp[27rif2jj(|/)] are the multiplicative periods, in the regime of y corresponding 


to small (pi,P 2 )- Plugging into (4.12) gives Fyac{^ 2 , 3 ) for arbitrary c; we then proceed to 
study this result at large c. 
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For two disjoint intervals and arbitrary n, the period matrix is 38 


, . n—1 


\ nj \ n J 


k=l 


(5.10) 


Specializing to n = 3, the period matrix is given by 


Q/ A ^ 2-Fi (3, 3; 1; 1 y) 

™ V3 2n(i.i;i;!/) 


(5.11) 


This is a highly symmetric genus-two Riemann surface: there is only a single modulus y, as 
opposed to the 3g — 3 = 3 moduli of a generic genus-two surface. 

To express the Schottky coordinates in terms of y, we need to invert the power series 


expansion given in (2.17). The fact that qn = ^22 implies that in Schottky coordinates, ^ 2,3 


has Pi = P 2 = p, as a quick inspection of (2.17) reveals. Our results are applicable when 


p 1, so (2.17) forces us to take qu 1 too. From (5.11), this is just the short interval 


limit, p <C 1, often studied in the context of 2D CFT Renyi entropy: taking p -C 1 in (5.11) 
yields multiplicative periods 


Qii 


qi2 






lOp^ 

7^ ^ 6 ^ 


+ 


27 

-15 

y 


2 p 


29y^ 

19683 

73%2 

’ 729 


+ 0{y^), 
41812/3 
6561 


(5.12) 


+ 0 (/) 


Finally, we obtain the series expansion of p and x in terms of y by inverting (2.17) using 
(5.12) and the explicit results for the coefficients c{n, m, |r|) and d{n, m, r) given in 40 . The 


result is 


v { y ) + 19683 


2/3 + 


26 


19683 


y^ + 


5768 


4782969 


27 


56 


28 


12892 


y 27^ 27^ 19683 


y - 


2/3 + 
3044 
6 ^ 


47429 

43046721 


2/3 + 


10582844 

10460353203 


y^ + 0(2/3). 


2/^+ 0(2/3), 

(5.13) 


Note that x{y) diverges linearly for small y 
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We can now compute the vacuum free energy F^g^dy) = — log .^vac(2/), and hence the 


We note that p{y) is nothing but the square of the larger eigenvalue of the Schottky generators themselves: 
eig[Li{y)) = p{yd‘^, where Li{y) = L 2 {y) are the Schottky generators in the y I regime. This was already 


computed in ^ 39 ,|^ ; see in particular equation (3.8) of [39| , with k = l,n = 3. One can then find 
x{y) using the Schottky relations. Such an algorithm is an alternative to that presented in the text. 


44 






























vacuum contributions to Renyi entropy, in a 1 short interval expansion, where 


OO 

ZUv) = E , 

hi,h2=0 


(5.14) 


We will further expand this result at large c and compare to those of 33 


In order to perform this expansion, we need to be a bit careful: because powers of x 


introduce inverse powers of y, it is not manifest in (5.14) that the short interval expansion 


can be meaningfully organized in powers of p{y). We need to know something about how 


Chi,h 2 {x) scales with large x, and hence small y. Fortunately, we can read this off from (4.25). 


Keeping terms to leading order in ?/ —)■ 0 at each order in 1/c, and ignoring coefficients, (4.25) 


and (5.13) imply that for h > hi > 0, 


lim lini p{y)^Ch^^h-hAy) ~ C>(/) + -0{y 

c-^oo y^O C 


h+2\ 




(5.15) 


vn=2 


Therefore, we can indeed ignore higher order terms in the sum over h when we expand in 
small y. 

Without further ado, the results are as follows. At £ = 1, 2 we find 


Fhac; l{y) 
Fhac; 2(2/) 


y 


+ 


56 y^ 


+ 


2189 2/6 


24668?/^ 85 

177147 ' 4782969 ' 129140163 ^ 1162261467 ^ ^ 

8y^ 8y’^ 11122 y^ 51818 y^ 


(5.16) 


+ 


+ 


+ 


387420489 129140163 94143178827 282429536481 


+ 0{y 


105 


Comparing to the results of [^|^, we find agreement through 0{y^).‘^‘^ 33 only computed 
through 0{y^), so our term at 0(jf) is new. 


At three-loop order, ( 5.4[ ) implies a nonzero result. Evaluating (5.4) for pi = P 2 = p{y) 
and X = x{y), we find 

iAac;3(2/) = y^^ 


As discussed around (2.15), the authors of 33 computed FVac ;3 through 0{y^) only, and 
found zero. We now see that the first contribution appears at It is remarkable 

that a computation through 0{y^^) using twist fields would not have revealed the nonzero 
result! This speaks to the different strengths of the twist field method and the sewing 
expansion that we have performed. 
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We can compare (5.16) directly to I 3 in 33 . The mutual information cf. (2.7), has an overall factor 


of 1/2 for n = 3; including the anti-holomorphic part, as they do in 33 , contributes an overall factor of 2, 
so the two factors cancel. 
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Finally, nonzero all-loop results at 0{y^‘^) follow from (5.5) and (5.6): 


^ _ 2/^^ /4(3704 + 590c + 1250^) 5^ 

c-j>z[y) - ^ 5c(5c + 22) } 


+ 0{y 


IS', 


(5.18) 


5.2.2 Other genus-two Renyi entropies 

Consider the three-interval Renyi entropy on the plane with n = 2. In this case, the replica 
manifold is a genus-two manifold characterized by three moduli that parameterize the 
positions of the intervals modulo conformal symmetry^ Our results for the free energy for 
P 11 P 2 ^ 1 and general x can therefore be regarded as (universal contributions to) Renyi 
entropies for the case of three disjoint intervals and n = 2. 

The period matrix of ^ 3^2 is known in terms of Lauricella functions 50 . To apply our 


results, one would first need to understand the relative spacings of intervals that corresponds 
to pi,p 2 1, by using the map from Schottky space to the period matrix. We do not pursue 
this geometric picture here. It is clear, however, that not all intervals need to be short, 
because x is allowed to be general. Thus, we have implicitly provided the first computations 
of universal contributions to 2D CFT Renyi entropies that do not require all intervals to be 
short. 

One can also consider the case of one interval on the torus with n = 2. The replica 
manifold has two moduli, namely, the temperature and interval length. Our methods can 
again be applied to this case to derive universal contributions to the Renyi entropy from 
the stress tensor sector. This has been done perturbatively in a high or low temperature 
expansion in using different methods that cannot access terms at two-loop and 

beyond in a large-c expansion, unlike the sewing method here. We note that pi and p 2 as 
a function of the moduli have been computed perturbatively in 15,51,52 . We leave the 
remaining explicit calculation for future work. 


5.3 The separating degeneration limit 


An important predecessor of the present work is 28 , where the relation between Zvac and 


3D gravity was first enunciated precisely. Yin tested this relation at genus two, focusing on 
the separating degeneration limit of the Riemann surface, where S becomes the union of two 
tori. Before we probe this region of moduli space with our new results, let us briefly review 


the work of 28 


^^Besides the case of n = 2 Renyi entropy for two intervals (for which the replica manifold is a torus with 
complex structure r given by a known function of the interval length [49| ), this is the only replica manifold 
that spans its entire genus g = (Y — l)(n — 1) moduli space. 
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For the sake of easy comparison to 28 , we use his notation in this subsection. We write 
the elements of the period matrix as 


ffii — p , ^22 — O' , r2i2 — y 


(5.19) 


We also dehne the multiplicative periods 


q = , s = , v = 2niy 


2'Kia 


(5.20) 


The separating degeneration limit corresponds to the limit u —)■ 0 with (g, s) hxed, where q 
and s parameterize the complex structure of the two tori. 

In the 1/c expansion, 


computed parts of FVac; 05 -^vac; i and Fvac ;2 at genus two using 


a variety of methods, all of which agree: 

• Demanding a match to the polar parts of extremal CFT partition functions at low val¬ 
ues of k = c/24, which are hxed by invariance under the genus-two modular group Sp(4, Z). 
That this match should hold follows from the dehnition of extremal CFTs, theories that 
have no non-trivial Virasoro primaries of dimension less than k + 1 above the vacuum. 


• For F/ac;i, direct computation of (5.7). 

• Direct computation of Zvac written as a sum over bilinears of torus one-point functions 
of Virasoro vacuum descendants. This is similar to what we do in the present work. 

Although it is not our focus here, Fq is given by a certain Liouville action whose origins 
we explain in Appendix]^ In order to write the expressions for F/ac; i and F^s.c; 2 , we must 
dehne the holomorphic Eisenstein series, normalized as 


EP = 




I _(jm 


(5.21) 


m=l 


For n = 2,4, these hatted versions relate to the usual Eisenstein series as 


Eo = -—^ q + 3g^ -|- 4g^ -|- 0{q^) 


K = 


24 

Ejjq) - 1 
240 


(5.22) 


q + 9q^+ 28q^+ 0{q^) . 


The results of 28 , which we denote are as follows: in the separating degeneration 
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limit f —?• 0, 


2q 


2s 


F^i = - E (riE^2 + Y3T-2 

n=2 ^ 

-^"*1 


- AE^E^ 


(5.23) 


y ^ - 2(g + s) + 45(g^ + s^) + 72qs + 745(g^s + qs^) + 3720g^s^j + 0(g^, /) 

0{v^) 


ail' 


c0 


^2 = 2 ^^ 


l-q 


-EP 


1 — s 


-E^, 


+ 2Av‘^ ( + s) - i^qs ) + 0(g'‘, s"^) ) + 0{v'^) 


13 1 


45 


36 8 


16 


(5.24) 


Note that these are non-perturbative in q and s throngh O(n^), and the leading term in 
is just the sum of one-loop free energies on two tori with periods q and s. (Note that 
in order to recover the (^(n*^) piece of one relies on (4.24).) Expanding everything 

through O(g^s^), we find 

^vlc; 1 = (q^ + q^ + S^ + S^) 

— (2qs{q -1- s -f 3gs)(2 -|- 3g -|- 3s -|- 8gs) j 

-f ^ — 2(g -I- s) -I- 45(g^ -f s^) -|- 72qs + 745(g^s -|- gs^) -|- 3624g^s^ j 

+ 0 (^;^g^s^) 


(5.25) 


and 


E:^, = 2Av^ [ q^s^ [l + Zq + Zs+^-qs 


2 2 2 


1 1 


3 2 

13 1 


45 


(5.26) 


+ 24n^ ( g^s^ ( ^ + ^(? + «) “ ) ) + C»(n^ g^ s^) 


^36 8'^ ^ 16 

We are now in a position to extend these results using our computations. As in the 


^^The semiclassical expansion of the free energy in 28 was performed in powers of 1/fc = 24/c. We expand 
in 1/c, and define according to the 1/c, rather than the 1/k, expansion. Thus, our equals 24 

times the S 2 found in 


28 
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previous subsection, our goal is to perform the map {pi,P 2 ,a;} h-)■ {q,s,v}, express F^ac in 
these variables, and expand at large c. 


In 28 , the following relations were established: 


P. = 9 (l - p"(24’) - P* + lmE2 - ^E; + f 

P2 = Pl(o- O p) . 

All we need now is to derive x{q,s,v). We do so by inverting one of the Schottky relations 


in equation (2.17), 


n+m 


= X 


+ x ^2 PiPT ^ d{n,m,r)x'^ . (5.28) 

n,m=l r=—n—m 

Note that for s = g = 0, we have x = . So we can write this as x = + 0(g, s, gn, sv). 

The hnal result for x is rather appealing. 


X 


= e^ - + v^) + 0{v^qs) . 


(5.29) 


We derive this in Appendix [C} 


Plugging equations (5.27) and (5.29) into (4.12) enables us to extend the results of 28 


m 


two ways. First, we can now give the (9(n^) part of the one- and two-loop free energies (5.25) 
and (5.26), respectively, through 0{q^s‘^,q'^s^), not only through O(g^s^). Second, and more 


importantly, we can write down some of the leading terms as n —)■ 0 for all loops. 
We hnd 


rp _ pYin 

-Tvac; 1 — -f^vac; 1 


+ ( qQS (g^(210 + 2764s + 11865s^) + s^(210 + 2764g + 11865g2)) + O(gS^) 


0 (n® 


(5.30) 


and 


^vac ;2 = EJ^^ “ (l4(g2 s^) 283gs(g -F s)) -h 0(g^s'^)) + 0{v^) . (5.31) 


To read off the free energy at three loops and beyond, we plug (5.29) into (5.4)-(5.6). At 
three-loop order, we hnd 


Fvac; 3 = q S 


44/133124 86528 6 


. + + O(n^) ) + 0(gS®^;^ g^sV) 

\ ZiO 10 


(5.32) 
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To derive terms at higher orders in v for hxed q and s, we would need to expand x beyond 
O(n^). Likewise, to derive terms at higher orders in q and s for hxed v, we would need to 
include more terms in the sewing expansion, like PiP 2 ^ 4 , 5 (a;). Finally, the all-loop expansion 
is completed by the term^ 




/4(3704 + 590c+ 125 c2) 


5c(22 + 5c) 


(5.33) 


6 Discussion 


We close with a discussion of some open questions, progressing from obvious directions for 
future work to the more speculative. Some directions for future work were mentioned in the 
text. 


In the realm of Renyi entropy, performing the calculation suggested in section 5.2.2 


for three intervals would give a satisfying derivation away from a short-interval expansion. In 
addition, one can straightforwardly apply our results to the case of the n = 2 Renyi entropy 
for a single interval on the torus, at least in a high or low temperature expansion. One 
need only perform the map between Schottky coordinates and the temperature and interval 
length; this map has already been partially performed in [l5|[5l|[^ . No results have yet 
been derived for the torus case beyond one loop. 


• One can consider including local operator insertions on S. The sewing procedure 
remains a sum over sphere correlation functions, now with these extra operator insertions. 
The operators generate non-vacuum states in the OFT. Taking S to be a replica manifold, 
one can thus compute excited-state Renyi entropies by the sewing procedure. Such entropies 
have been computed in CFT using twist-field and holographic methods (e.g. 18,19,53-^). 
As we have tried to demonstrate, the sewing construction is likely to provide a complementary 
approach that operates at hnite c, so this seems like an especially worthwhile pursuit. It 
would be easy, for instance, to read off the 0{c^) terms from the above procedure: these 
would be predictions for bulk one-loop corrections to the Einstein action evaluated on the 
“punctured handlebody”. 


18,19,53 


• It would be nice to prove that nowhere in the moduli space, except at the separat¬ 
ing degeneration point, does the genus-two partition function truncate in a 1/c expansion 
(whereas our method could only access the regime of small pi,P 2 )- This seems highly likely 
to be the case. Understanding the structure of the Schottky sum rules in Appendix [C] could 


^^This result should be contrasted with footnote 6 of 


28 
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also be enlightening. 


• In our analytic bootstrap of section we could equally have used Virasoro conformal 
blocks rather than global conformal blocks. In this case, the crossed blocks are related to 


the original blocks by the fusion and braiding matrices. These are known in closed form 56 


so our conclusions can also be phrased in terms of OPE coefficients of Virasoro primaries 
rather than quasi-primaries. The Virasoro approach is in principle more efficient, as it will 
§x the four-point function in terms of even fewer pieces of data, and it would be worthwhile 
to make this precise. An interesting demonstration of this fact comes by way of the W 3 


correlator (yVWWW) ^ as computed in (3.41)-(3.42): up to the norm of IT, the Virasoro 
approach would hx iyVWWW) without having to compute even a single OPE coefficient. 

• We briefly considered CFTs with higher-spin symmetry; it would be straightforward 
to extend our computation of Zvac to higher orders for such theories. A more exciting 
prospect would be to compute the partition function on S in the presence of insertions 
that carry higher-spin charge. There is natural motivation for this from holography. In 
particular, while much work has been done to construct solutions of higher-spin gravity with 


nonzero higher-spin charge and solid-torus topology 57 , there has been no work on building 


solutions of higher-spin gravity of higher genus and with nonzero higher-spin helds turned 
on. A subset of such “higher-spin handlebodies” would be saddle points of the Euclidean 
higher-spin gravitational path integral with replica boundary conditions and nonzero higher- 


spin charge 48 ; accordingly, their action would be expected to match GET computations of 


Renyi entropy in states with higher-spin charge and/or chemical potentials. This calculation 
would be analogous to the one peformed in in the spin-2 case. Constructing such Renyi 
entropies via partition functions on replica manifolds endowed with higher-spin charge, rather 


than via twist helds 58 or Wilson lines 21 , would be an interesting application of the replica 


trick to the higher-spin setting. One might also try to make contact with the “spin-3 entropy’" 


of 59 . 


• Our results can be used to test the idea that Liouville theory provides an effective 
description of irrational CFTs with large central charge (see e.g. for a recent rehnement 
of this idea and references to earlier work). In particular, the 1/c expansion of the genus-two 
partition function can be checked against a diagrammatic calculation in Liouville theory. 


Upon hrst glance, the relation between 


CXD 


■^vac; 1 — / , 

hi,h2=0 


p{‘p 5* lim 


( 6 . 1 ) 


and the bulk graviton determinant (5.7) seems opaque. Nevertheless, these two quantities 
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are equal. Both formulae are written in terms of Schottky data, so it should be possible to 
hud a clean mapping between them. This would be a useful stepping stone to writing down 
a closed formula for the two-loop contribution to the bulk partition function, in analogy to 
the determinant (5.7). In the sewing prescription, the two-loop result simply requires us to 
sum over the 0{l/c) parts of Chi,h 2 {^) instead of just the (!7(c°) parts. Is there an equally 
simple prescription in the bulk, and if so do the primitive elements of the Schottky group 
play a privileged role as they do at one loop? 


• Part of the motivation for the present work was the one-loop exactness of the pure- 
gravity partition function on the solid torus. The current understanding of this result relies 
on an elegant and simple argument about Virasoro representation theory, which can be 
understood holographically. It can be derived without recourse to CFT by computing the 
energies of bulk excitations, or equivalently, by quantizing the phase space given by two 
copies of diff S'^/S'L(2, M) [^. Still, it would be very satisfying to derive this result from a 
more direct perspective in the bulk. For example, while the solid torus partition function of 
a pure higher-spin theory is also believed to be one-loop exact, we do not know the analog 
of diffS*^ in that context; there should be a more direct argument one can make in the bulk. 
Understanding this exactness from the perspective of the bulk diagrammatic expansion could 
provide insights useful for higher genus. 

On the other hand, a perhaps cleverer approach would be to derive the partition function 
from the S'L(2,M) x SL{2,M.) Chern-Simons formulation of 3D gravity. Einstein-Hilbert 
gravity and Chern-Simons theory are non-perturbatively inequivalent, but it is believed that 
the semiclassical expansion around a well-dehned saddle point can be performed in either 
formulation. The Chern-Simons approach builds in the topological nature of 3D gravity, 
whereas the loop expansion of 3D gravity in the metric formulation is no simpler than it is in 
higher dimensions, despite the absence of propagating bulk degrees of freedom. Presumably, 
such a computation would be manifestly one-loop exact, in analogy to similar truncations in 
compact Chern-Simons theory [^. (More precisely, all higher-loop effects could be absorbed 
in a renormalization of the Newton constant.) A Chern-Simons approach would also have 
the beneht of immediately generalizing to pure higher-spin gravity. The challenge to carrying 
this out is that both the gauge group and the topology are non-compact. There has been 
progress in recent years in computing Chern-Simons partition functions for non-compact 
gauge groups (see e.g. [^), but the requisite technology does not yet exist for the solid 
torus. This technology would represent a signihcant advance in our understanding of 3D 
gravity. 


52 



Acknowledgments 


We thank Matthias Gaberdiel, Jared Kaplan, Christoph Keller, Albion Lawrence, David 
Poland, David Simmons-Duffin, Herman Verlinde, Roberto Volpato and Xi Yin for helpful 
discussions. MH, AM, and EP wish to thank the Aspen Center for Physics for hospitality 
during this work, which was supported in part by National Science Foundation Crant No. 
PHYS-1066293. MH and ICZ were supported in part by the National Science Foundation 
under CAREER Crant No. PHYlO-53842. EP was supported in part by funding from the 
European Research Council under the European Union’s Seventh Framework Programme 
(FP7/2007-2013), ERC Crant agreement STC 279943, “Strongly Coupled Systems”, and 
in part by the Department of Energy under Crant No. DE-FC02-91ER40671. ICZ was 
supported in part by the Department of Energy under Crant No. DE-SC0009987. AM is 
supported by the National Science and Engineering Research Council of Canada. 


A An OPE limit of Chi,h2{^) 

The goal in this Appendix is to show that the x —)■ 0 behavior of Ch^^h-hiix) is given, as in 


(4.25), by 


hm \imCh,,h-hA^) ~ 0{x-^) + -0(x-"+2) + ( V 4 1 
c^oo c \ • 

\n=2 


(A.l) 


where we again have written h 2 = h — hi, and we restrict to h > hi > 0. We are ignoring hi- 
and h-dependent coefficients at each order, and displaying only the leading singular behavior 
at each order in 1/c. 


The upshot is that (A.l) follows from considering the t-channel OPE limit of the four- 


point functions that dehne the Chi,h 2 {x)- The expansion (A.l) follows from the c-scaling 
of the OPE coefficients and norms that appear in the conformal block decomposition. The 


leading 0{c^) term in (A.l) arises from identity exchange; the leading 0{l/c) term, from T 
exchange; and all higher order terms in 1/c, from exchange of all other quasi-primaries in 
the Virasoro identity representation. 

Let us give more detail. Recall that the Chi,h 2 {x) are dehned in terms of sums over 
four-point functions of vertex operators. 


where the Hilbert subspaces Hh, are spanned by operators of holomorphic dimensions hj. 
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Vertex operators and V*”(0, z) are just chiral CFT operators 'ip{z) and 

respectively, dressed with 2 ;-dependent factors. In the 2 ; —)■ 0 limit, the dressing factors are 


hnite (cf. section 4.1), so we can ignore them. Thus, the x —)■ 0 limit of Ch^^h 2 {.x) is simply 
the t-channel limit of a weighted sum over four-point functions of CFT operators, including 
descendants. 

A four-point function of pairwise identical quasi-primary operators of dimensions hi and 
h 2 can be written in a global conformal block decomposition as 


O',; 


(V’(oo)0(x)0(l)^/>(O)) = X ^ V] -^x'^° 2 Fi{ho, ho] 2ho] x 


(A.3) 


where C^^o are OPE coefficients, and h = hi + ^ 2 . We are expanding this correlator in 
the X —)■ 0 channel. Four-point functions involving secondary operators can be written 
using derivatives acting on an expression of the above form. For our purposes, we only 
allow Virasoro descendants of the identity to run in the internal channel: (P G {1, T, A,...}. 


Indeed, for the purposes of establishing (A.l), consideration of the exchange of these three 


operators alone will be sufficient: that is, we associate the scaling in (A.l) with specific terms 


in (A.3). Let us write out the first three terms in (A.3) coming from the Virasoro identity 


block, O e {1,T,A}: 
(V’(oo)0(x)0(l)V'(O)) = 


X 


— h 


/^2 _ 


xSEi(2,2;4;x) + a^SEi(4,4; 8; x) + 0(x®)^ 


(A.4) 


We have substituted the explicit operator norms. 

When -0 and 0 are in the same global conformal family, fusion onto the identity is allowed, 
and C^^i 7 ^ 0 and independent of c. This yields a term of order O(c^) and 0(x“^). For 
a given {hi,h 2 ), there is always such a term in the definition of Ch^^h^ix), since the latter 
are defined as a sum over all correlators involving operators at levels (hi, h 2 )- This is most 
obvious when hi = h 2 , because Chi,h 2 {x) will include correlators of four identical operators; 
but even when hi 7 ^ h 2 , the definition of Chi,h 2 {x) includes correlators of arbitrary derivatives 
of T, all of which are in the same global conformal family. For instance, 6 * 2 , 4 (a^) includes 
(T(oo) d‘^T{x) d‘^T{l) T(0)), which permits fusion onto the identity when x —)• 0. There are 


generically no cancellations among terms in the sum (A.2). This accounts for the first term 


on the right-hand side of (A.l). This leading behavior was also observed in |4b) . 

The second term in (|A.l ), at (P(l/c), comes from the second term in (A.4). Because 


is c-independentp^ and N't = c/ 2 , this term contributes at 0 {l/c) compared to the 


26 


For instance, when tp = (p is a. quasi-primary, = /lA//. 
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identity exchange, but not beyond. As explained above, the dehnition of Chi,/i 2 (a;) always 
includes such terms. This accounts for the second term on the right-hand side of (A.l). 

The hnal terms in (A.l), at OiXjc^') and beyond, come from exchange of the level-four 
quasi-primary A in ( |A.4 ). Because its inverse norm has an inhnite expansion in 1/c, this 
will contribute a term to all orders in a 1 /c expansion, thus accounting for the 

remaining terms in (A.l). 


B More details on the sewing construction 

B.l Operators and norms 

In this section, we list the operators and their norms at the hrst six levels of the Virasoro 
vacuum representation. To ensure that we have not missed any, it is useful to expand the 
holomorphic Virasoro vacuum character, Xvac^ 


Xvac = Trvac(g^‘’ 


= g-^/24 


n 

n=2 



_ ^-c/24^^ + ^2 ^ ^3 ^ + 2q^ + ...) . 


(B.l) 


One can branch Xva.c into global S'L( 2 ,]R) characters, thereby counting the number of quasi¬ 
primary helds. The resulting generating function, call it Xqp, is 

Xqp = - l){l - q) ^ + 2q^ + .... (B.2) 


In terms of the degeneracy d{h) of all level h operators, the degeneracy of level h quasi¬ 
primaries is d{h) — d{h — 1) (for h > 1). We use the shorthand O = O{0) to denote 
operators. 

• Level 2: there is one quasi-primary operator, the stress-energy tensor, T = L_2|0), 
with norm At = c/2. 

• Level 3: there is one secondary operator, = dT = L_3|0), with norm Aog = 2c. 

• Level 4: there are two operators, 

L_4^ |0) , of) = L_4|0) . (B.3) 
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The operator is the commonly studied quasi-primary often denoted A, and 
is secondary. Their norms are 


( 2 ) 




A/T( 2 ) = 5c . 


o 


(B.4) 


Level 5: there are two operators , 


O™ = L_1 I - ii-4 110). Ot' = i-5|0) , 


)( 2 ) _ 


(B.5) 


where both of them are secondary. Their norms are 


M^(i) - 4 c ( c + y j , 


A/1 ( 2 ) = 10c . 


o 


(B.6) 


Level 6: there are four operators (all acting on |0)), 


- 1L_4L_2 + 1L_3L_3, 


20 


63 


9 


9 


(B.7) 


( 2 ) (60c+ 78) 3(42c + 67) 93 , r r r 


(70c+ 29)' 


(70c + 29) 


(70c+ 29)' 


Oq — L_iL_i (L_ 2L_2 — - L_4) , 
^6 — -^- 6 ) 


where and are quasi-primary, and (9® and Og"^^ are secondary. Their norms 


are 


.r 4 3 (2c-1) (5c+ 22) (7c+ 68) 

= hh ^ + 29) > = - c- 


A 63 

^r 4,0 f 22\ 35 

A/1{3) — 72 c ( c + — , A/1(4) — —c . 


(70c + 29) 


5/’ 2 


The secondary operators L_„|0), n > 2 have the form ‘^^T/{n — 2)!, with norms n(n^ — 
l)c/12. 


B.2 Four-point functions of vertex operators 

In this section we provide more details on the transformation properties of the vertex opera¬ 
tors which were used in computation of the four-point functions Chi,h 2 {x) in subsection 
The final expressions for Chi,h 2 {x) in terms of x are reported in the main text and are not 


4.2 
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repeated here. 

We will need the explicit expressions for the vertex operators at inhnity. For h = 2,3, 
we have 


W“*(T,cx)) = lim z^T{z) 

Z^OO 

W“*(C> 3 , oo) = lim f - z^dT{z) - 4z^T{z)] . 

For h = 4, we have 

W“*(A,oo) = lim z^A{z) 

Z^OO 

V°^\oP,oo)= lim (lz^d^T{z) + 5z'^dT{z) + 10z^T{z)) . 

For h = 5, we have 

oo) = lim f - - z^L^) 0^^\z) 

Z^OO \ J 

oo) = lim f - -z^^d^T{z) - 3z^d^T{z) - l3z^dT{z) - 2{)z'^T{z)] . 
z^oo V 6 / 

Finally, for h = 6, we have quasi-primary vertex operators 

«d) = lim of(z), t = {1,2} , 

Z^OQ 

and secondary vertex operators 

(/“‘(Of ,oo) = lim [z^^ + z'^L, + —L'f) of'(z), 

Z^QO \ Z / 

V°^\0^^\oo) = lim {^z^^d*T{z) + ]:Z^^d^T{z)+ z^^ d‘^T{z) 

z^oo \24 0 Z 

+ 33z^ dT{z) +33 z^T{z)^ . 

With these in hand, we start with Ch, 2 {,x). Using the dehnitions 

V°^\T,x) = {x-lfT{x) 

W”(T,1) = (x- 1)2T(1) 

W«(C),0) = 0(0) 


(B.8) 


(B.9) 


(B.IO) 


(B.ll) 


(B.12) 


(B.13) 
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that follow from section 4.1, Ch, 2 {x) takes the general form 

d{h) 


cux) = - 1)‘(v'”“(o‘‘',« d) t(x) r(i) oi‘>( 0 )). 

i=l '* ' 

We next consider the functions Ch, 2 ,(x)- For h = 3 we have 


V°^\0,,x) W"(03,l) >, (B.15) 


where cxd) is given in (B.8) and 


l/°“'(C>3,a;) = -{x - If dT{x) - A{x - If T{x), 
W”(C> 3 ,1) = (x - l)=^aT(l) - 4 (a: - 1)2T(1) . 


(B.14) 


(B.16) 


The expressions for Cx^^fx) and C^^fx) can then be easily obtained using (B.16) and the 
vertex operators given above. 


C Schottky parameters in the separating degeneration 
limit 

Here, we provide details of the map between Schottky space and the period matrix in the 


separating degeneration limit considered in section 5.3 


{Pi,P2,a:} {^,5,-?;}. 

The hnal result, perturbative in v but non-perturbative in q and s, is 

Pi = 9 ("l - - V* 

P2 = Pl(o- -f-)- p) 

x = e^ - + vf + Oiy^qs) . 


(C.l) 


(C.2) 


The hatted Eisenstein series were dehned in (5.21). The hrst two relations were derived 
; here, we will derive the last one. 


m 


The computation is an exercise in series solutions of algebraic equations. We make a 
series ansatz for x, 

OO 

j=o 


(C.3) 








plug this and the expressions for pi and p 2 in (C.2) into the perturbative Schottky relation 


oo n+m 

= X + X ^ PiP^ ^ d{n,m,r)x'^ , (C.4) 

n,m=l r=—n—m 

and solve order-by-order for Xj{q,s). 

An immediate question that may occur to the reader is how we are able to obtain a result 
that is non-perturbative in q and s, despite only having access to d{n,m,r) to hnite 
order in {n,m). The answer is that we were able to infer various sum rules obeyed by the 
d{n,m,r) that we believe to hold for all {n,m): 



Ml 

o 

(C.5a) 

r=—n—m 

n+m 

d{n, m,r) r = 0 

(C.5b) 

r=—n—m 

n+m 

d{n, m, r) = 0 

(C.5c) 

r=—n—m 

oo n+m 

(C.5d) 

n,m=l r=—n—m 

n+m 

X", d{n,m,r)+ = 0 

(C.5e) 

r=—n—m 

where E 2 was dehned in (5.21). We have also found a set of sum 

rules obeyed by the 

c{n,m, r ) that appear in the other two Schottky relations (2.17): 

n+m 

XI c(n, m, r) = 0 , {n,m) ^ (0,0) 

r=—n—m 

n+m 

XI c(n, m, r ) = 0 , n 7 ^ 0 

(C. 6 ) 

r=—n—m 

00 m 

X S™ X = ^^2 • 

m=l r=—m 

We have checked all of these identities through m = n = 7 using 

the tables of 40 ^ 

Actually, we have proven (C. 6 ), as well as (C.5a)-(C.5c). Proof of (C. 6 ) follows from compar- 


ing a series solution for pi and p 2 using (2.17) to the known solution (C.2), and demanding 


are grateful to the authors of 40 for sharing the relevant Mathematica notebooks. 
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consistency throngh O(n^). Proof of (C.5a)-(C.5c) follows from demanding that all three 


pertnrbative relations in (2.17) yield the same resnlt for {pi,p 2 ,a^}: hence, having proven 


(C.6), we can nse these to derive snm rnles obeyed by d{n,m,r). Proof of (C.5d) and (C.5e) 


is nndonbtedly possible nsing similar methods. 


With these snm rnles in hand, we proceed to invert (C.4). At 0{v^), we mnst solve 


n+m 


1 = Xo{q,s) + Xoiq,s) ^ ^ d{n,m,r)xo{q, sY 

n,m=l 


(C.7) 


r=—n—m 


Bnt (C.5a) implies that xo{q^s) = 1 solves (C.7). At 0(n), we mnst solve 


n+m 


1 = a;i(g, s) + a;i(g, s) d{n,m,r)r 

n,m=l r=—n—m 


(C.8) 


This time, (C.5b) implies that the second term vanishes, leaving xi{q^s) = 1. The analysis 
at O(n^) is nearly identical, and (C.5c) implies X 2 {q,s) = 1/2. 


At O(n^), the hrst non-trivial snm appears in the series expansion: 


1 1 

^ = x3(g,s) + ^ 

n,m=l 


n+m 

q"‘s'^ ^ d{n,m,rY 

r=—n—m 


(C.9) 


Plngging in (C.5d) leads to 


xsiq, ^ - 4 ^ 2^2 


Finally, at we mnst solve 


(C.IO) 


- oo ... oo 

^4(g,s) + - g”s™^d(n,m,r)r3 + - Y 


4! 


n,m=l 


n,m=l 


(C.ll) 


The hnal snm rnle (C.5e) eliminates the last term, and (C.IO) leaves ns with 


Xi{q, s) = ^ - 4^2^2 


Pntting this all together, we hnd the advertised resnlt in (C.2). 


(C.12) 


We believe that the above snm rnles may have interesting applications in other stndies 
of genus-two Riemann surfaces. It would be interesting to understand, for instance, why the 


sum (C.5d) factorizes. A systematic exploration of all sum rules obeyed by these coefficients 


would be worthwhile. It seems likely that higher order sum rules may be expressible in terms 


of holomorphic Eisenstein series (5.21) 
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D The order-c part of the free energy and the sewing 
construction 


In subsections 2.2| and 4.1, we reviewed the sewing construction, which expresses the partition 
function of an arbitrary CFT on a genus-^f Riemann surface in terms of 25 f-point functions 
on the sphere, as illustrated in hgure However, the formulas in those subsections, such as 
(4.4), only give the order-c° and higher (in 1/c) terms in the free energy F = — InZ; they 


miss the order-c term. (This term depends on the full metric on the Riemann surface, not 
just its complex structure; in other words it depends on the choice of representative of the 
Weyl class.) This is adequate for the purposes of this paper, since our main interest is in 
the higher-order terms in 1/c. However, for completeness, in this appendix we will explain 
how to obtain the order-c term within the context of the sewing construction. 

As an illustrative example, consider an arbitrary CFT on a flat torus with modular 
parameter r. The partition function is well-known to be 


Z{t) = {pp) = (pp) ^ d(h, h)p' 


V 


(D.l) 


h,h 


where p := and d{h, h) is the multiplicity of operators of weights h, h. To calculate 


Z{t) using (4.4) (more precisely, its generalization including the antiholomorphic sector). 


we need to compute the coefficient Applying the dehnitions (4.5) and (4.20), we hnd 
simply Cf^j^ = d{h, h). Hence (4.4) would give Z{t) = ^ d{h, h)p^p^] thus, we are missing 

the factor (pp)“'^/^^. 


We proceed with a brief recap of the sewing construction]^ For convenience we will 
assume that the metric on our Riemann surface M is smooth. We now cut it along a circle 
and glue in two disks, which we call Di and D 2 , to obtain a new manifold Mq (which may 
be connected or disconnected). We choose the metric on these disks in such a way that the 
metric on the new surface is still smooth. This implies that Di and D 2 can be glued together 
to make a sphere with a smooth metric, which we’ll call S. We consider coordinates Zi ^2 on 
Di 2 which, when extended to S, obey zi = \lz 2 - 

The path integral on M can be computed by inserting a complete set of states on the 
circle where it has been cut. By the state-operator mapping, this is equivalent to inserting 
a complete set of operators on D 12 at 2 : 1^2 = 0, with an appropriate inverse metric (to 


^®We will closely follow the discussion of the sewing construction in sections 9.3 and 9.4 of 62 . However, 


that reference considered CFTs with vanishing total central charge (in the context of string theory), so the 
issue we are focusing on here did not arise there. When comparing to that reference, note also that the 
correlators there are unnormalized, whereas (as throughout this paper) ours are normalized. 
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be determined below) on the space of operators: 


Z(M) = Z(Mo) , 


(D.2) 




where the superscripts denote that Oi is inserted at the origin of the zi and Oj at the origin 
of the Z 2 coordinate system. More generally, we can start with arbitrary operators Oa - ■ ■ 
on M\ 

Z{M){Oa ■ ■ ■)m = Z{Mo) G'HOa ■ ■ ■ . (D.3) 




To £x the inverse metric we consider the case where M happens to include a patch 
D[ that is diffeomorphic to Di, with an operator Ok inserted at = 0. Cutting M along 
the boundary of D[ and gluing in Di and D 2 yields Mq = M U S, where the S is covered by 
coordinates z'^ and Z 2 = Equation (D.3) then becomes: 




Z{M)Z{S) Y, G'HOa ■ ■ ■ 0Y)M{0Y0iYs . (D.4) 


Mo 
{^2) 




For this to hold for arbitrary Ok and arbitrary insertions Oa - ■ ■, it must be that 

G,i = Z(S){0‘"‘>0p>)s = Z(S)5.j , 


(D.5) 


where Qij is the Zamolodchikov metric. 


Now that we have hxed G®-^, the partition function (D.2) becomes 


z(M) = ^)(^5;e«(op'op>)„.. 


Z(S) 




(D.6) 


It is often useful to add a parameter p to the sewing construction, so that the coordinate 
identihcation is Z 1 Z 2 = p- (Even though p can be absorbed in a coordinate transformation 
on Mq, it is useful to £x the coordinate system on Mq and use p to vary the modulus of M.) 
This can be done by replacing Z 2 by Z 2 in the formulas above, and dehning Z 2 = p^ 2 - 
have O^Y = pGpGoY, so 

Z(M) = . (D,7) 


Cutting M along g non-contractible cycles, where g is its genus, reduces it to a sphere. 
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This yields the formula (4.4), except with a product of sphere partition functions Z{Si) ■ ■ ■ Z{Sg) 
in the denominator. The free energy on any sphere is proportional to c, so these factors con¬ 
tribute such a term to F{M). (The coordinate transformation from local coordinates zi^2 in 
the vicinity of each operator insertion to the single coordinate 2 ; covering the plane leads to 
the definition of the operators V°^^, Id®"' explained in subsection 4.1 ) 


To illustrate the application of (D.7), let us return to the example of the flat torus. Set 
{3 = Imr, and let the horizontal cycle have circumference 27r; thus the total area is We 

will cut it along the horizontal cycle. For Di and D 2 we use unit round hemispheres. Thus S' 
is a round unit sphere, while Mq is a cylinder of circumference 27r and length 27r/3 with round 
endcaps. In the next paragraph we will compute the ratio Z(Mq)/Z(S) using the Liouville 
action, finding precisely the prefactor appearing in the expression (D.l) 

for the torus partition function. 

In order to compute Z(Mo)/Z(S'), we will compute the change in the partition function 
Z{Mq) under a small change in f3, and then integrate the result up from (3 = 0 (noting 
that Mo 1/3=0 = S). Under a Weyl transformation, the partition function gets 

transformed by the Liouville action: 


Z = e^^Z, 


Sr = 


2471 


(g'^’^daUidbU + Ru 


(D.8) 


We will let ds^ be the metric with cylinder length 27i(3, and ds^ with cylinder length 27i{f3 + 
6(3). Hence the Weyl transformation relating them is close to the identity, with u of order 
6(3, and we can work to first order in uj. Since the cylinders have the same circumference, uj 
vanishes on the cylinder, uj can also be taken to vanish on, say, the bottom endcap, while on 
the top endcap it transforms the hemisphere into a hemisphere attached to a thin cylinder 
of height 2776 (3. On this endcap, R = 2, so 



j ^/§RuJ = J v^2cn = J Va - J \/I = 47r^d/3 . 

(D,9) 

Hence 

Sl = '^6(3. 
b 

(D.IO) 

Integrating from (3 = 

0 , we hnd 



TTC 

\nZ{Mo) = \nZ{S) + —^, 

6 

(D.ll) 


as promised. 
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